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Nomenclature 


a 


d<P> 


L 

m12 

"FI 


vector  of  pitch  line 

major  and  minor  axes  of  gear  cradle  ellipse 

major  and  minor  axes  of  pinion  cradle  ellipse 

gear  head  cutter  diameter 
Mean  pitch  cone  distance 
gear  ratio:  m^  “ 

pinion  cutting  ratio 


gear  cutting  ratio 


n 


(F) 


q 

qF*  qp 


unit  normal  to  pinion  and  gear  tooth  surface 
direction  angle  of  normal  in  coordinate  system 
parameters  of  pinion  and  gear  machine-tool  setting 


(P)  (F) 

r'  ,  r 
c  c 


radius  of  generating  surface  measured  in  plane  X 

(i  -1,2) 


(i) 


m 


=  0 


V  up 


(F)  (P) 

V  9  V 


(1)  (2) 

v  ,  v 


~A 


v 

~c 


generating  cone  surface  coordinate 
surface  E  ,  E  contact  point  velocity 

r  Jr 

surface  E^,  contact  point  velocity 

velocity  of  intersection  point  at  gear  cutter  axis  in  plane  II 
velocity  of  intersection  point  at  pinion  cutter  axis  in  plane  II 


gear  cutter  width 
mean  spiral  angle 

pinion  and  gear  pitch  angle 

pinion  root  angle 

orientation  angle  of  pinion  and  gear  cradle  ellipses 


v 


£ 


V 


rotation  angle  of  frame  relative  to  frame  about  axis 

rotation  angle  of  frame  S,  relative  to  frame  S 

f  n 

pinion  and  gear  dedendum  angle 
pinion  machine-tool  settings 

generating  pinion  and  gear  cone  surface  coordinates 

motion  parameter  of  cradle  ellipse 

motion  parameter  of  pinion  and  gear  ellipses 


0) 


(2) 


plane  of  normals 

pinion  and  gear  tooth  surfaces 

generating  surfaces  of  pinion  and  gear 
generating  surfaces  rotation  angle 

pinion  and  gear  blade  angle 
pinion  and  gear  angular  velocities 

cradle  angular  velocities  for  cutting  the  pinion  and  gear 


Cartesian  Coordinate  Frame 

connected  to  tool  cone,  j  =  F,  p 
s 

S(J)  conntected  to  cradle 

c 

connected  to  machine  frame,  i  =  1, 


m 

Sj,  connected  to  pinion,  gear 

S,  fixed  to  machine,  used  for  mesh  of 


2 


I-  and  E, 
F  1 


vi 


sf  fixed  to  frame  of  gearbox,  used  for  mesh  of  E  and  E 

1  P  2 

Sn  connected  to  plane  of  normals 


SUMMARY 


Spiral  bevel  gears  are  used  in  many  applications  where  mechanical 
power  must  be  transmitted  between  intersecting  axes  of  drive  shafts. 

Namely,  two  such  applications  are  the  rear  axle  differential  gearbox  for 
land  vehicles  and  the  transmissions  used  in  helicopters.  For  spiral 
bevel  gears,  there  is  a  continuing  need  for  ever  stronger,  lighter  weight, 
longer-lived  and  quieter  running  gears.  Above  all,  a  rapid  and  economical 
manufacturing  method  is  essential  to  the  industries  that  use  bevel  gearing 
in  their  products. 

For  many  years,  the  Gleason  Works  (Dudley,  1962;  Anon.  Gleason  Works, 
1964  and  1980)  has  provided  the  machinery  for  manufacture  of  spiral  bevel 
gears.  There  are  several  important  advantages  to  the  Gleason  methods  of 
manufacture  over  hobbing  methods.  The  machines  are  rigid  and  produce  gears 
of  high  quality  and  consistency.  The  cutting  methods  may  be  used  for  both 

milling  and  grinding.  Grinding  is  especially  Important  for  producing 

/ 

hardened  high  quality  aircraft  gears.  Both  milling  and  grinding  are 
possible  with  Gleason's  method.  The  velocity  of  the  cutting  wheel  does  not 
have  to  be  related  in  any  way  with  the  machine's  generating  motion. 

Generally  speaking  Gleason's  method  for  generation  of  spiral  bevel 
gears  does  not  provide  conjugate  gear  tooth  surfaces.  This  means  that  the 
gear  ratio  is  not  constant  during  the  tooth  engagement  cycle,  and, 
therefore,  there  are  kinematical  errors  in  the  transformation  of  rotation 
from  the  driving  gear  to  the  driven  gear.  The  research  that  had  been 
performed  by  the  Gleason  Works  was  directed  at  the  minimization  of  gear 
kinematical  errors  and  the  improvement  of  gear  bearing  contact  by  using 
special  machine-tool  settings.  The  determination  of  such  machine-tool 
setting  is  accomplished  by  a  computer  program.  It  is  known  that  Oerlicon 
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(Switzerland)  and  Klingelnberg  (West  Germany)  have  developed  methods  for 
generation  of  spiral  bevel  gears  that  can  provide  conjugate  gear  tooth 
surfaces.  The  disadvantage  of  these  methods  is  that  the  gear  tooth 
surfaces  cannot  be  ground  and  the  tooth  element  proportions  are  unfavorable 
due  to  the  constant  height  of  the  teeth. 

The  objective  of  the  new  method  for  generation  of  spiral  bevel  gears 
presented  herein  was  to  find  a  way  to  eliminate  the  kinematical  errors  for 
spiral  bevel  gears,  obtain  gears  with  higher  contact  ratio,  and  improve 
bearing  contact  and  conditions  of  lubrication,  while  using  the  existing 
Gleason's  equipment  and  retaining  all  the  advantages  of  the  Gleason  system 
of  manufacturing  such  gears.  The  proposed  method  for  generation  is  based 
on  the  following: 

(i)  Four  surfaces  -  two  generating  cone  surfaces  (E  and  E_)  and  gear 

P  F 

and  pinion  tooth  surfaces  (^2  and  E^ )  are  in  continuous  tangency  at 
every  instant.  The  ratio  of  angular  velocities  in  motion  of  the 
above  mentioned  surfaces  satisfies  the  requirement  that  the 
generated  pinion  and  gear  transform  rotation  with  zero  kinematical 
errors . 

(ii)  The  cones  have  a  common  normal  at  the  instantaneous  point  of  contact 
but  their  surfaces  interfere  with  each  other  in  the  neighborhood  of 
contact  point. 

(iii)  The  point  of  contact  of  the  above  mentioned  surfaces  moves  in  a 

plane  (II)  that  is  rigidly  connected  to  the  gear  housing.  The  normal 
to  the  contacting  surfaces  lies  in  plane  II  and  performs  a  parallel 
motion  in  the  process  of  meshing. 

(iv)  Due  to  the  elasticity  of  gear  tooth  surfaces  their  contact  is 

spread  over  an  elliptical  area.  The  proposed  method  for  generation 
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provides  that  the  instantaneous  contact  ellipse  moves  along  but  not 
across  the  gear  tooth  surface.  This  provides  improved  conditions 
for  lubrication  and  a  higher  contact  ratio  will  result. 


(v)  Until  now,  the  reduction  of  kinematical  errors  of  Gleason  spiral 
bevel  gears  was  a  subject  of  the  computer  search  for  the  optimal 
machine-tool  settings.  The  computer  program  developed  for  this 
purpose  is  called  the  TCA  (Tooth  Contact  Analysis)  program.  The 
proposed  method  is  for  direct  determination  of  machine-tool  settings 
that  result  in  zero  kinematical  errors  because  the  gear  tooth 
surfaces  are  generated  as  conjugate  gear  tooth  surfaces. 

(vi)  A  new  TCA  program  directed  at  the  simulation  of  bearing  contact  and 
the  influence  of  errors  of  assembly  and  manufacturing  has  been 
developed. 

The  contents  of  this  report  covers  the  new  method  for  generation  of 
spiral  bevel  gears,  their  geometry,  the  bearing  contact  and  simulation  of 
meshing.  Special  attention  has  been  paid  to  the  proposed  principle  of 
performance  of  parallel  motion  by  two  related  ellipses  that  results  in  the 
desired  parallel  motion  for  the  contact  normal. 

1.  THE  GLEASON  MANUFACTURING  METHOD 

The  gear  cutter  cuts  a  single  space  during  a  single  index  cycle.  The 
gear  cutter  is  mounted  to  the  cradle  of  the  cutting  machine.  The  machine 
cradle  with  the  cutter  may  be  imagined  as  a  crown  gear  that  meshes  with  the 
gear  being  cut.  The  cradle  with  the  mounted  head  cutter  rotates  slowly 
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about  its  axis,  as  does  the  gear  vhich  is  being  cut.  The  combined  process 
generates  the  gear  tooth  surface.  The  cradle  only  rotates  far  enough  so 
that  one  space  is  cut  out  and  then  it  rapidly  reverses  while  the  workpiece 
is  withdrawn  from  the  cutter  and  indexed  ahead  in  preparation  for  the 
cutting  the  next  tooth.  The  desired  cutter  velocity  is  provided  while  the 
cutter  spins  about  its  axis  which  itself  moves  in  a  circular  path. 

We  consider  that  two  generating  surfaces,  Ip  and  Ip,  are  used  for  the 
generation  of  the  pinion  tooth  surface,  Ej ,  and  the  gear  tooth  surface, 
respectively.  Both  sides  of  the  gear  tooth  are  cut  simultaneously  (duplex 
method)  but  both  sides  of  the  pinion  tooth  are  cut  separately  (single 
method).  The  basic  machine-tool  settings  provide  that  four  surfaces,  Ep, 
Ip,  Ij  and  I2,  are  in  contact  at  the  main  contact  point.  In  the  process  of 
meshing,  surfaces  Ip  and  1^,  and  respectively  surfaces  Ip  and  ^  contact 
each  other  at  every  instant  at  a  line  (contact  line)  which  is  a  spatial 
curve.  The  shape  of  the  contact  line  and  its  location  on  the  contacting 
surfaces  is  changed  in  the  process  of  meshing.  The  generated  pinion  and 
gear  tooth  surfaces  are  in  contact  at  a  point  (contact  point)  at  every 
instant . 

A  head-cutter  used  for  the  gear  generation  is  shown  in  Fig.  1.1.  The 
shapes  of  the  blades  of  the  head-cutter  are  straight  lines  which  generate  a 
cone  while  the  head-cutter  rotates  about  axis  C-C.  The  angular  velocity 
about  axis  C-C  does  not  depend  on  the  generation  motion  but  only  on  the 
desired  cutting  velocity.  Two  head-cutters  are  used  for  the  pinion 
generation;  they  are  provided  with  one-sided  blades  and  cut  the  respective 
tooth  sides  separately.  The  head  cutter  is  mounted  to  the  cradle  of  the 
machine.  Fig.  1.2  shows  schematically  the  positioning  of  the  cradle  of  the 
gear  machine,  the  head  cutter  and  the  gear  to-be  generated. 


2.  GENERATING  SURFACES  AND  COORDINATE  SYSTEMS 


The  generating  surface  is  a  cone  surface  (Fig.  2.1).  This  surface  is 

generated  in  coordinate  system  while  the  blades  of  the  head-cutter 

c 

rotate  about  axis  C-C  (Fig.  1.1).  The  generation  of  gear  tooth  surfaces  is 

based  on  application  of  two  tool  surfaces,  I_  and  Z_,  which  generate  gears 

r  Jr 

1  and  2,  respectively.  The  generating  surfaces  (generating  cones)  do  not 

/ p\  C  P) 

coincide:  they  have  different  cone  angles  ipc  and  ^  *  and  different 

(  F  )  (  P  ) 

mean  radii  r';  and  r^  (Fig.  2.1, a).  Special  machine-tool  settings,  AE. 

C  1 

and  Alx  (Fig.  2.3,b),  must  be  used  for  the  generation  of  the  pinion. 
Considering  the  generation  of  gear  2  tooth  surface  we  use  the 

(p) 

following  coordinate  systems:  (i)  S  which  is  rigidly  connected  to  the 

(2) 

generating  surface  Z  (Fig.  2.1,b);  (ii)  the  fixed  coordinate  system  S 

p  m 

whicn  is  rigidly  connected  to  the  frame  of  the  cutting  machine,  and  (iii) 

the  coordinate  system  S2  which  is  rigidly  connected  to  gear  2  (Fig.  2.2). 

(2) 

In  the  process  of  generation  the  generating  surface  rotates  about  the  X 

m 

(p) 

-  axis  with  the  angular  velocity  w  while  the  gear  blank  rotates  about 

(2)  (2) 

the  Z2  axis  with  the  angular  velocity  w'  .  Axes  X^  and  intersect 

each  other  and  form  the  angle  90°  +  ^2  ~  ^2'  w*iere  ^2  *s  t*ie  dedendum  angle 

(2) 

for  gear  2.  Axis  Xm  is  perpendicular  to  the  generatrix  of  the  root  cone 
of  gear  2.  The  coordinate  system  Sf  shown  in  Fig.  2.2  is  rigidly  connected 
to  the  housing  of  the  gears  and  will  be  used  for  the  analysis  of  conditions 
of  meshing  of  the  gears. 

Considering  the  generation  of  the  pinion  we  use  the  following 

(F) 

coordinate  systems:  (i)  Sc  which  is  rigidly  connected  to  the  generating 
surface  (ii)  ^  which  is  rigidly  connected  to  the  frame  of  the 
cutting  machine  and  (iii)  which  is  rigidly  connected  to  the  pinion  (gear 
(Fig.  2.3).  Axes  and  Z^  do  not  intersect  but  cross  each  other;  AE^ 


and  AL^  are  the  corrections  of  machine-tool  settings  which  are  used  for  the 

improvement  of  meshing  of  the  gears.  In  the  process  of  generation  the 

generating  surface  rotates  about  the  -axis  with  the  angular  velocity 

(F) 

a)  while  the  gear  1  blank  rotates  about  the  Z^  -axis  with  the  angular 
velocity  W(1).  Axes  and  form  the  angle  90°  -  t  where  A  is 

the  dedendum  angle  of  gear  1;  axis  is  perpendicular  to  the  generatrix 

of  the  root  cone  of  gear  1. 


3.  GENERATING  TOOL  SURFACES 

The  tool  surface  is  a  cone  and  is  represented  in  the  coordinate  system 
as  follows  (Fig.  2.1) 


r^^cotif/^  -  u  cos/^ 

C  C  5  yc 

u  .sinij/^  ^sin0, 
j  c  j 

u.sin<l/^  ^cos0. 

J  c  J 


(j  -  F,  P) 


(3.1) 


where  Uj  and  0^  are  the  surface  coordinates. 

The  coordinate  system  (j  =  F,  P)  is  an  auxiliary  coordinate 

system  which  is  also  rigidly  connected  to  the  tool  (Fig.  2.1,b).  To 
represent  the  generating  surfaces  and  Ep  in  coordinate  system  ^ 
use  the  following  matrix  equation 


x(J)_ 

c 

s 

y(J) 

C 

-  [M(j)] 
cs 

y(J> 

s 

z(^ 

c 

s 

_  1 

1 

1 

0 

0 

0 

x(j)' 

s 

0 

cosq^ 

+  sinq . 

J 

+b^sinq^ 

y(j) 

s 

0 

+sinqj 

cosqj 

b .  cosq , 

3  j 

2<j> 

s 

0 

0 

0 

1 

1 

Here:  Ik  and  q^  are  parameters  which  determine  the  location  of  the  tool  in 

the  coordinate  system  \  Henceforth,  the  upper  sign  corresponds  to  the 
generation  of  a  left-hand  spiral  bevel  gear  that  is  shown  in  Fig.  2.1  and 
lower  sign  for  a  right-hand  spiral  bevel  gear. 

Equations  (3.1)  and  (3.2)  yield 


C3>  -  r(3>cot*l)>  -  u.coS„c(3) 


c 

,(3> 

c 

,<J>  = 


)t^ 
c  c 


u.sintjr  J 'sin(  0.+  q,  )  +  b.sinq 
j  c  j  J  j  J 

u  .sim//^cos(  0  ,+  q  .  )  +  b  cosq 
j  c  J  J  J  J 


(3.3) 


where  j  =  (F,  P). 

The  unit  normal  to  the  generating  surface  Zj  (j  =  F,  P)  is  represented 
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by 


n(J> 

~c 


N(j) 

~c 


I  Sc 


3r«>  3rP> 

7-rr—  ,  where  —  x  - - - 

(j)|  ~c  30.  3u . 


(3.4) 


J 


,(j) 


Using  Eqs .  (3.3)  and  (3.4)  (provided  u.sin^J  ^  0),  we  obtain 

1  c 

=  sinij/^i^^  +  cos^^[sin(0.+q  )j^^  +  cos(  0.  -  q.)k^^] 
~c  c  ~c  c  3  j  ~c  J  J  ~c 


(3.5) 


4.  EQUATIONS  OF  MESHING  BY  CUTTING 

Generation  of  .  We  derive  the  equation  of  meshing  of  the  generating 

surface  ^  and  the  gear  tooth  surface  using  the  following  procedure: 

Step  1:  First,  we  derive  the  family  of  surfaces  ^  that  we  represent  in 

the  coordinate  system  Such  a  family  is  generated  while  the 

coordinate  system  S^F^  is  rotated  about  the  -axis  (Fig.  2.3).  We 

c  m 


recall  that  the  generating  surface  £_  is  rigidly  connected  to  S 

*  ( 

coordinate  transformation  in 
the  following  matrix  equation 


(F) 


The 


coordinate  transformation  in  transition  from  S^F^  to  is  represented  by 

c  m 


X(1) 

x(F) 

m 

c 

m 

m  c 

,<F) 

c 

,<l> 

Z(F) 

m 

c 

_i  . 

. 

(4.1) 


Here  (Fig.  2.3): 
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1  0 


[mOKP)]  -  o  cos<j>  sin<}>  0 

me  r  r 


(4.2) 


0  -sintL,  cost})-  0 


0  0 


where  <f>  is  the  angle  of  rotation  about  the  Xv  -axis. 
"  m 

Using  Eqs.  (4.1),  (4.2)  and  (3.3),  we  obtain 


r<F)cot*<F>  -  Vo8^F> 


Ir<»]  - 


upsin^  ;sinTp  +  bpsin(qp  +  4^) 


(4.3) 


where 


T  =■  0  r  q  _  +  <t> 

F  F  +  hF  F 


(  r  I 

UpSiniJ;  cosxp  +  bpcos(qp  ^  <j)p) 


The  upper  sign  corresponds  to  the  right-hand  spiral  bevel  pinion. 

Equations  (4.3),  with  parameter  4>p  fixed,  represent  a  single  surface  of  the 
family  of  generating  surfaces. 

Step  2:  The  unit  normal  to  the  generating  surface  Ip  may  be  represented  in 
the  coordinate  systems  as  follows: 


,  where  Nv 


(4.4) 


We  may  also  use  an  alternative  method  for  the  derivation  of  the  unit 
normal.  This  method  is  based  on  the  matrix  equation. 


[n(1)]  =  [L(1)(F)]  [n(F)] 


(4.5) 


Matrix  [l/*^F^]  may  be  determined  by  deleting  the  4^  column  and  rows  in 
me' 

matrix  (4.2). 

(F) 

The  column  matrix  [n  ]  is  given  by  vector  equation  (3.5). 

c 

After  transformations,  we  obtain 


[n(l>] 

m 


.  (F)  . 

cosi^c  sintp 


(4.6) 


.(F) 

COSlp  COST— 
c  * 


Step  3:  We  derive  the  equations  of  the  relative  velocity,  v  ,  as 

- i - 

follows : 


V<F1> .  V<F> .  vo> 

~m  ~m  -m 


(4.7) 


where  v^F^  is  the  velocity  of  a  point  N  on  surface  and  ^  is  the 


velocity  of  the  same  point  N  on  surface 
(F) 

Vector  v  is  represented  by  the  equation 
~m 


v(F>  ,  JF>*r<» 

~m  ~m  ~m 


(4.8) 


Here  (Fig.  2.3): 


I 

■ph  mi  i  mi 
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- : -  - - - - - : : 1; 1  *  *•  **  *  ■  * * !  *  ** 


(4.9) 


-to 


(F) 


Vector  r^^  is  represented  by  equation  (4.3). 
Equations  (4.8)  and  (4.9)  yield 


r  (F ) i  =  (F) 

lv  J  O) 
m 


.(1) 


m 


-y. 


(i) 


m 


(4.10) 


Gear  1  rotates  about  the  Z^  -  axis  with  the  angular  velocity  w 
(Fig.2.3).  Since  u/ ^  does  not  pass  through  the  origin  of  the 


(1) 


(1) 


we  substitute 


(1) 


m 


U) 


coordinate  system  S 

through  0(1)  and  the  vector  moment  represented  by 


by  an  equal  vector  which  passes 


m 


n( 1 )n  (1) 
0  0.  x  0) 

m  n  ^ 


Then,  we  represent  v^^  as  follows 

~m 


m  n  ~m 


(4.11) 


-m 


-m  ~m 


Here: 
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[v<mP,)l 


-(j/1)(y^1)  +  AE1)cos(y1  -  Aj) 

U)(F)zO)  +  w(1){(x(l)  _  LsinA1)co8(Yl  -  Aj) 
-w(F)yi1)  +  a)(1)(y^U  +  AE1)sin(Y1  -  Aj) 


-  (z^°  +  AL1)sin(Y1  -  Aj)} 


(4.14) 


Step  4;  The  equation  of  meshing  by  cutting  Is  represented  by 


„<»  .  v<F1>  -  0 
~m 


(4.15) 


Using  equations  (4.15),  (4.14),  (4.6)  and  (4.3),  we  obtain 

j-uF  +  [r^F^cot^F^  -  LsinAj  -  ALj  tan(yj  -  Aj)]cos^F^|  sinTp  + 

b  [siniJ;(F)sin(q  +  <j>  )  +  cos^<F)sin0  ^ - 

F  c  F+  F  -  rc  F  cosCYj  -  Ap  J 


AEj^sint^1^  -  cosif/F\an(Y  j  -  A^coslp] 


fl(uF’  0F*  V  "  0 


(4.16) 


“fi  3  "ny 

0) 


This  equation  relates  the  generating  surface  coordinates  (up,  9p)  with  the 
angle  of  rotation  (4>p). 


iV  Vv*’  V.* 


Generation  of  ®y  using  a  similar  procedure  we  may  obtain  the  equation 

of  meshing  for  surface  Zp  and  surface  as  follows 

Step  1:  Equations  of  the  family  of  generating  surfaces  £p  represented  in 

(2) 

the  coordinate  system  S  are: 

m 


fx(2)l 

m 

r^P^COtlj/P^  -  UpCOS||/P^ 

m 

3 

UpSin^p)sinTp  +  bpsin(qp  ;  <J>p) 

z(2) 

m 

^UpSin^p)c°sTp  +  bpcos(qp  +  4>p) 

(4.17) 


(2) 

Step  2:  The  unit  normal  to  £  is  represented  in  S  by  the  column  matrix 
- * -  p  m 


(P) 


sin^ 


cos>J^P^  sinT 

c  t 


cos4^P^  COST 

c  p 


(4.18) 


where 


T  »  0  +  q  +  0 

P  P  P  P 


The  upper  sign  corresponds  to  the  left-hand  spiral  bevel  gear  and  the  lower 

to  the  right-hand  spiral  bevel  gear. 

(P)  (2)  (P2)  (2) 

Step  3:  The  velocities  v  ,  v  and  v  are  represented  in  S  as 
- -  ~ra  ~m  ~m  m 


follows 


Iv«p>, 

m 


(4.19) 


(2) 

ym  cos(y2  -  A2) 

(2) 

■(x''  ’  +  LsinA0)cos(y0 
m  z  1  z 


y<2)S1„(Y2  -  42) 


zl2>  sln(Y2  -  V 


-  h£p>>  -  tv<2>: 


Step  4:  The  equation  of  meshing 


n<2)  .  v(P2)  =  0 
~m  ~m 


yields  the  following  equation 


(4.20) 


(4.21) 


f2(up»  Qp*  4>p)  3 

[up  -  (r^P^cotij/ +  LslnA2)cos^P^  ]sinTp  + 


(P)  (P)  mP2  ~  ® ir'(Y2—  ) 

bpsln(qp  +  Vsln*c  +  bPC0S*c  sln0p  —cos(7^-a2) - 0 


(4.22) 


Here: 


P2  (2) 

oj 


S.  ORIENTATION  OF  THE  PINION  CRADLE 

Henceforth,  we  will  consider  two  auxiliary  coordinate  systems  and 
S^,  that  are  rigidly  connected  to  the  gear  and  pinion  cutting  machines, 
respectively  (Fig.  5.1). 

Figure  2.2  shows  coordinate  system  that  is  rigidly  connected  to  the 

(2) 

gear  cutting  machine  and  to  coordinate  system  Sm  .  Coordinate  system 

is  also  rigidly  connected  to  the  housing  of  the  gear  train  and  the  meshing 

of  the  generated  pinion  and  gear  will  be  also  considered  in  coordinate 

system  S^.  Axis  is  the  instantaneous  axis  of  rotation  of  the  pinion  and 

the  gear  -  the  pitch  line  (the  line  of  tangency  of  the  pinion  and  gear 

pitch  cones).  The  origin  0^  of  coordinate  system  coincides  with  the 

(2^ 

points  of  intersection  of  3  axes:  ’ ,  Z ^  (Fig.  2.2)  and  z^  (Fig.  2.3). 

(2) 

Here:  X  is  the  axis  of  rotation  of  the  gear  cradle;  z~  is  the  axis  of 
m  Z 

rotation  of  the  gear  being  in  mesh  with  the  generating  gear  Zp  and  Z^ 
is  the  axis  of  rotation  of  the  pinion  being  in  mesh  with  the  generating 
gear  Zp  and  the  gear  member. 

Figure  2.3  shows  coordinate  system  that  is  rigidly  connected  to  the 

pinion  cutting  machine  and  coordinate  system  The  origin  0,  of 

m  n 

coordinate  system  coincides  with  the  origin  0^  but  the  orientation  of 
system  with  respect  to  represents  a  parameter  of  the  machine-tool 
settings  (proposed  by  Litvin,  1968).  This  parameter  is  designated  by  € 
in  Fig.  5.1. 

The  coordinate  transformation  in  transition  from  to  is  based  on 
the  following  considerations. 

(i)  Consider  two  auxiliary  coordinate  systems,  and  ,  that  are 

rigidly  connected  to  systems  and  S^,  respectively  (Fig.  5.1, a  and 
Fig.  5.1,b),  axes  Z.  and  z,  coincide  with  the  pinion  axis.  Initially 
coordinate  system  coincides  with  and  system  with  S^. 
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(il)  Assume  now  that  coordinate  systems  and  are  rotated  about  the 
Z  -axis,  the  axis  of  the  pinion,  through  the  angle  £  (Fig.  5.1,c). 

This  angle  determines  the  orientation  of  coordinate  system  with 
respect  to  S^,  or  the  orientation  of  the  pinion  cutting  machine  with 
respect  to  the  gear  cutting  machine.  Matrix  [M^]  represents  the 
coordinate  transformation  in  transition  from  to  and  is  represented  by 
the  following  equation 

'"fh1  ’  '"fa1  lMab'  'V 

r  2  .  ,  2 

cosecos  Yj  +  sin  Yj 
=  -sinecosYj 

cosy^siny^  (1  -  cose) 

(5.1) 


sinecosy^  cosy^siny^l  -  cos£)  0 

cose  sinesiny^  0 

2  2 

-sinesiny^  cosesin  y^  +  cos  y^  0 

0  0  1 


6.  PLANE  OF  NORMALS 

It  will  be  proven  below  that  the  generating  surfaces  Zp  and  Zp  contact 
each  other  at  a  point  that  moves  in  the  same  plane,  II.  The  generated 
pinion  and  gear  surfaces,  Z^  and  Z^,  also  contact  each  other  at  every 
instant  at  a  point  that  coincides  with  the  point  of  tangency  of  surfaces  Z^ 
and  Zp.  However,  we  have  to  emphasize  that  surfaces  Zp  and  Z^ 
(respectively,  Zp  and  Z^)  are  in  line  contact  and  the  instantaneous  line  of 
contact  moves  over  the  contacting  surfaces.  The  special  property  of  the 
method  developed  for  generation  of  spiral  bevel  gears  is  that  the  contact 
point  moves  in  a  plane  (plane  II)  that  is  rigidly  connected  to  the  fixed 


coordinate  system  S^.  Also,  the  common  normal  to  the  contacting  surfaces 
lies  in  plane  n  and,  as  it  will  be  shown  later,  performs  a  parallel  motion 
in  the  process  of  meshing.  Plane  II  is  called  the  plane  of  normals  and  it 


is  determined  as  the  plane  that  passes  through  the  instantaneous  axis  of 
rotation,  Z^,  and  the  normal  N  to  the  generating  surface  Zp  (Fig.  6.1). 

The  above  normal  passes  through  point  N  that  is  the  main  gear  contact 
point . 

Figure  6.2  shows  the  orientation  of  coordinate  system  with  respect 

to  coordinate  system  S^.  Origin  0^  coincides  with  origin  0p  and  axis 

coincides  with  the  Z  -  axis.  The  coordinate  transformation  in  transition 

n 

from  to  is  given  by  the  matrix 


cosn 

-sinn 

0 

tw  - 

sinn 

cosn 

0 

(6.1) 

0 

0 

1 

The  determination  of  angle  n  is  based  on  the  following  considerations: 
(1)  unit  vector  i  can  be  represented  in  coordinate  system  by  the 
following  matrix  equation 


cosr)  -sinn 

sinri  cosn 

0  0 


0 

0 

1 


1 

0 

0 


cosn 

sinn 

_0 


(6.2) 
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(2)  The  unit  normal  n^  to  the  generating  surface  Ep  is  represented  as 
follows 


>4p>>  -  ^<2)"»:2,i 


(6.3) 


Here  (Fig.  2.2) 


r  (2)  , 

tLfm  1 


COSA2  0  -sinA2 


1  0 


sinA^  0  COSA2 


(6.4) 


(2) 

1 

and  (4.18)  yield 


Unit  vector  n  is  given  by  column  matrix  (4.18).  Equations  (6.3),  (6.4) 


(p)  (p) 

sin^c  cosA2  -  cos^c  cosTpsinA2 


(P) 

cos^  sinT 
c  P 


(P)  (P) 

sin^c  8inA2  +  cosipc  cosTpcosA2 


(6.5) 


(3)  Vectors  i^n^  and  n^P^  are  mutual  perpendicular,  i.e.,  i£n^  *  nf^  =  0. 


This  yields  that 


tanh 


(P) 

tan:|r  COSA2  -  coslpsinA2 
sini„ 


(6.6) 


Here: 


6p  +  <Jp  +  K 


with  the  upper  sign  for  the  left-hand  gear, 

(P) 

=  otp  for  the  convex  gear  tooth  side, 

(p) 

=  180°  -  otp  for  the  concave  gear  tooth  side 
where  otp  is  the  angle  of  the  cutter  blade. 

Equation  (6.6)  determines  the  angle  of  orientation  n  of  the  plane  of 
normals  (Fig.  6.2). 

7.  PERFORMANCE  OF  PARALLEL  MOTION  OF  A  STRAIGHT  LINE  PROVIDED  BY  TWO 
RELATED  ELLIPSES. 

An  important  part  of  the  proposed  approach  is  a  new  technique  directed 
at  the  performance  of  a  parallel  motion  of  a  straight  line  provided  by  two 
related  ellipses.  By  using  this  technique  it  becomes  possible  to  provide  a 
parallel  motion  for  the  common  normal  to  the  gear-pinion  tooth  surfaces. 
This  motion  is  performed  in  a  plane  (the  plane  of  normals)  that  is  rigidly 
connected  to  the  gear  housing  and  has  the  prescribed  orientation. 

It  is  well  known  that  a  translational  motion  of  a  straight  line  may  be 
performed  by  a  parallelogram  linkage  (Fig.  7.1, a).  Consider  that  a 
straight  line  slides  by  its  points  A  and  C  along  two  circles  of  equal 
radii.  Vectors  v^  and  v^  which  represent  the  velocities  of  points  A  and  C 
of  the  moving  straight  line  are  equal.  The  moving  straight  line  AC  being 
initially  installed  parallel  to  the  center  distance  OD  will  keep  its 
original  direction  in  the  process  of  motion. 

The  discussed  principle  of  translational  motion  of  a  straight  line  may 
be  extended  for  the  case  where  the  straight  line  slides  along  two  mating 
ellipses  (Fig.  7.1,b).  These  ellipses  have  the  same  dimensions  and 
orientation  and  again  the  velocity  vectors  v^  and  v^,  are  equal.  Consider 
that  the  moving  straight  line  is  initially  installed  parallel  to  OD  where  0 
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AV, 


W,V 


and  D  are  the  foci  of  symmetry  of  the  ellipses.  Then,  with  v^  =  the 
moving  straight  line  will  keep  its  original  direction  in  the  process  of 
motion. 

Figure  7.1(a)  and  Fig.  7.1(b)  show  a  translational  motion  of  a  segment 
of  a  straight  line  (AC)  with  constant  length.  A  more  general  case  is 
represented  in  Fig.  7.1(c).  The  straight  line  slides  over  two  ellipses 
whose  dimensions  and  orientation  are  different.  The  length  of  segment  AC 
which  slides  along  the  ellipses  is  changed  in  the  process  of  motion.  The 
problem  is  how  to  provide  a  parallel  motion  of  the  the  moving  straight 
line.  We  call  this  motion  a  parallel  one  because  the  straight  line  has  to 
keep  its  initial  parallel  to  line  OD  where  0  and  D  are  the  foci  of  symmetry 
of  two  mating  ellipses.  Unlike  the  cases  which  are  shown  in  Fig.  7.1(a) 
and  Fig.  7.1(b)  the  motion  of  straight  line  is  not  translation  because  the 
velocities  yA  and  v^,  of  the  tracing  points  A  and  C  are  not  equal.  The 
distance  between  the  sliding  points  A  and  C  is  changed  in  the  process  of 
motion.  It  will  be  proven  that  the  required  parallel  motion  of  straight 
line  may  be  performed  with  certain  relations  between  the  dimensions  and 
orientation  parameters  of  two  guiding  ellipses  which  are  shown  in  Fig. 
7.1(c). 

Consider  that  an  ellipse  (Fig.  7.2)  is  represented  in  coordinate 
system  by  the  equations 


Xa  =  °»  ya  =  aPC0SV  Za  =  bP  sinyP 


(7.1) 


where  ap  and  bp  are  the  lengths  of  the  semimajor  and  semiminor  axes,  respec¬ 
tively,  D  is  the  origin  of  coordinate  and  the  symmetry  focus  of  the  ellipse; 
parameter  up  determines  the  location  of  a  point  on  the  ellipse. 


m 

I 


Coordinate  system  has  the  same  origin  as  Sa  and  the  orientation  of 
with  respect  to  S&  is  given  by  angle  6p.  To  represent  the  ellipse  in 
coordinate  system  we  use  the  following  matrix  equation. 


[r(P)]  =  [M  ][r<P)] 

r»  tl3.  cl 


(7.2) 


where 


1  0 


[M  ]  =  0  cos6p  -sinfip  0 


0  sinSp  cos6p  0 


0  0 


(7.3) 


Matrix  equation  (7.2)  yields 


/ p\  (  P) 

xn  =  °*  yn  =  apcos6pcosup  -  bpsin6psinyp 


=  apsin<5pcosyp  +  bpcos6psinyp 


(7.4) 


The  symmetry  focus  of  the  mating  ellipse  is  point  0  given  by 

coordinates:  C.  ■  0,  C„  =  y^°\  C,  =  z^  (C„  and  C.  are  algebraic 
i  i  n  j  n  i.  i 

values).  The  orientation  of  the  ellipse  is  given  by  angle  <$p  (Fig.  7.2). 
Equations  of  the  mating  ellipse  are  represented  in  coordinate  system  as 
follows 


*<F)  -  0,  ,<F> 

n  ’  Jn 


aFcos6FcosyF  -  bFsin6FsinyF  +  , 


ZR  =  aFsin6FcosyF  +  bFcos6FsinyF  +  C3 


(7.5) 


Consider  point  A  of  ellipse  1  is  determined  by  parameter  y 


Point  C 


PO* 

of  ellipse  2  is  determined  with  parameter  ypQ.  Henceforth  we  will  consider 
such  a  motion  where  Pp  -  Up^  =  Up  -  PpQ  =  y.  Line  AC  is  drawn  through 
points  A  and  C.  Our  goal  is  to  provide  that  line  AC  will  be  parallel  to 
the  center  distance  DO  for  any  value  of  the  motion  parameter  y  and  perform 
a  parallel  motion.  The  above-mentioned  goal  can  be  achieved  with  certain 
relations  between  the  ellipse  parameters  ap,  bp,  6p,  ap,  bp  and  6p.  We 
start  the  derivation  of  these  relations  by  considering  the  vector  equation 

DA(y)  +  AC(y)  =  D0(y)  +oc(y)  (7.6) 


Vector  equation  AC(y)  =  ADO  is  satisfied  with  any  value  of  y  if  the 

parallel  motion  of  line  AC  is  provided.  (Where  A  is  the  constant  required 

to  make  the  two  vectors  of  equal  magnitude.) 

Equation  (7.6)  yields 


(OC  -  DA)  •  jn 

(OC  -  DA)  •  k 
~n 


(A-l) (DO  ■  jn) 

(A-1)(D0  •  k  ) 
~n 


cos  q 
sin  q 


(7.7) 


Here:  i  and  k  are  the  unit  vectors  of  axes  Y  and  Z  ,  q  is  the  angle 

in  -n  n  n 

formed  by  axis  and  vector  DO  (Fig.  7.2).  For  the  further 

transformation  we  will  represent  y  and  y  as  follows 

P  F 


Up  =  Upo  +  U, 


Equations  (7.7), 


(b^sinq  - 


y  =  y  +  y 
mF  mF0  m 


(7.4),  (7.5) 
ausinq  -  b2l 


and  (7.8)  yield 
cosq  +  a2^cosq)cosy 


(7.8) 


+  (-b^siHq  +  b22cc>sq  +  ausinq  -  a22cosq)siny  =  0 


(7.9) 
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Here: 


all  =  apcos^pcosppo  -  bpsin5psinypo 

al2  =  apcos<Spsin^JP0  +  bpsin6pcosUp0 

a21  =  apsin6pcosV‘po  +  bpcos5pslnVpo 

a22  =  apsin6psin^po  “  bpcos<SpcosUp0  (7.10) 


Coefficients  b^,  b12,  t>21  and  b22  have  similar  expressions. 

Equations  (7.9)  must  be  satisfied  for  any  value  of  the  motion  parameter  y. 
This  means  that  equation  (7.9)  will  be  satisfied  for  any  value  of  y  if  the 
following  two  equations  are  satisfied  simultaneously. 


Hsinq  -  b21cosq 

-  a^sinq  +  a^cosq  *  0 

(7.11) 

i^sinq  -  b22cosq 

-  a^sinq  +  a22cosq  =  0 

(7.12) 

Equations  (7.10),  (7.11)  and  (7.12)  yield  a  system  of  two  pseudo-linear 
equations  in  two  unknowns  (cosy^  and  siny^): 

apsin(q-(5F)cosyF0  -  bFcos(q-6F)sinyFQ  =  dj  (7.13) 

bFcos(q-<5F)cosJiF0  +  aFsin(q-<5F)sinyF0  =  d2  (7.14) 

Here: 

d{  =  apsin(q-6p)coslip0  -  bpcos(q-6p)sinVJp0  (7.13) 

d2  =  apSin(q-6p)sinUp0  +  bpcos(q-6p)cosyp0  (7.16) 
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The  solution  of  equations  (7.13)  and  (7.14)  for  cosu^  and  siny  is 

ro  fo 

as  follows 


cosypo 

sinUpo 


aFdlsln^q-<^  +  bFd2c0S^“^p^ 

apSin2(q-6F)  +  b2cos2(q  -  6p) 

aFd2Sin^q"^F^  “  bpd1cos(q-6F) 
a2sin2(q-6F)  +  b2cos2(q-<5F) 


Equations  (7.13)  and  (7.14)  yield 


(7.17) 


(7.18) 


a2sin2(q-6F)  +  b2cos2(q-6p)  =  d2  +  d2 


(7.19) 


Equations  (7.19)  relates  3  parameters  of  the  mating  ellipse:  a„,  b 

r  F 

and  6  (Parameters  ap,  bp,  5p  of  the  first  ellipse  and  q  are  considered 
as  given).  Thus  Eq.  (7.19)  can  be  satisfied  with  various  combinations  of 
ap,  bp  and  6p. 

aV 

Example:  a  =  2,  b  =  1.25,  6  =  20°,  q  =  270°,  —  ■  3,  y  =  yFQ  *  0. 

F 

Then  we  obtain:  aF  =  2.2753,  <S F  =  34.3113°  .  These  ellipses  are  shown  in 
Fig.  7.2. 


A  linkage  which  may  perform  the  described  parallel  motion  of  line  AC  is 

shown  in  Fig.  7.3.  Link  1  is  in  contact  with  two  perpendicular  slots  and 

DN  =  a  ,  DM  =  b  where  2a  and  2b  are  the  lengths  of  the  axes  of  ellipse 
P  P  P  P 

which  is  traced  out  by  point  A.  Similarly,  link  1'  is  in  contact  with  two 
other  perpendicular  slots,  ON’  =  a_  and  0M'  =  b„,  point  C  traces  out 
ellipse  2.  The  orientation  of  a  pair  of  slots  with  respect  to  the  other 
one  is  determined  by  the  angle  (6„  -  6  )  (see  Fig.  7.2).  Links  1  and  1' 

f  p 

rotate  with  the  same  angular  velocity  to  =  Line  AC  will  perform  a 

parallel  motion  if  the  ellipse  parameters  satisfy  equation  (7.19). 
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8.  GEAR  MACHINE-TOOL  SETTINGS 


Input  Data.  The  following  data  are  considered  as  given  to  set  up  the  gear 
machine-tool  settings. 


N2,  the  teeth  number  of  the  gear 
Y2»  the  gear  pitch  angle  (Fig.  2.2) 

A2»  the  gear  dedendum  angle  (Fig.  2.2) 

Bp.  the  mean  spiral  angle  (Fig.  2.1) 

L  ,  0^Mo,  the  mean  distance  of  the  pitch  cone  (Fig.  2.2) 

W  ,  the  point  width  of  the  gear  cutter  (Fig.  8.1) 
the  blade  angle  of  the  gear  cutter  (Fig.  8.1) 

Gear  cutting  ratio.  The  gear  cutting  ratio  represents  the  ratio  between  the 
angular  velocities  of  the  cradle  and  the  generated  gear  and  is  designated 
by 


m  „ 
P2 


Jp) 


(8.1) 


(2) 

Henceforth,  we  will  consider  the  two  coordinate  systems,  S  and  Sf, 

m  i 

shown  in  Figure  2.2.  Axis  is  the  pitch  line  of  the  mating  spiral  bevel 
gears.  It  is  also  the  instantaneous  axis  of  rotation  of  the  pinion  and  the 
gear  that  transform  rotation  with  constant  angular  velocity  ratio.  Assuming 
that  the  generating  cone  surface,  Zp,  the  gear  tooth  surface,  Z2»  anc* 
pinion  tooth  surface,  Z^,  are  in  continuous  tangency  at  every  instant,  it 
is  required  that  the  instantaneous  axis  of  rotation  by  the  gear  cutting  must 
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coincide  with  the  pitch  line,  axis  Z 


Thus 


-  >>  .  „(2)  .  ^ 


This  means  that  vectors  ’  and  are  collinear. 


Vectors  u/P\  u/2  ^  and  k^  are  represented  in  coordinate  system  S 


as  follows 


«Jp)-  0 


(2)  (2) 

OJ  =  OJ 


'-sin(y2  -  A 2)' 


cos(y2  -  A  )_ 


0 


Equations  (8.2)  and  (8.3)  yield  that 


-uj(p)  +  u/2^sin(Y2  -  A2) 
-u(2)cos(y2  -  A2) 


Equation  (8.4)  results  in  that 


Main  Gear  Contact  Point.  The  main  gear  contact  point  is  the  center  of  the 


bearing  contact  on  the  gear  tooth  surface.  The  location  of  this  point  can 


be  represented  by  the  surface  coordinates  of  the  tool  cone  6p  and  up. 


We  can  vary  the  value  of  up  by  keeping  0p  fixed  and  obtain  the  desired 


location  of  the  main  gear  contact  point. 


We  start  with  the  case  when  Mq  coincides  with  the  pitch  point  Mq 


(Figure  8.2).  In  this  case  N  on  the  cutter  surface  is  the  contact  point 


of  the  pinion  and  gear  tooth  surfaces  since  the  normal  at  this  point 


intersects  the  instantaneous  axis  of  rotation,  Z^.  The  basic  relations 


for  the  left-hand  gear  for  this  case  are  as  follows: 


9  =90°  -  6  +  q  =T  +q 

P  P  P  P  P 


(8.6) 


LcosA2Sin(0p  -  qp)  LcosA^cosB^ 


(8.7) 


cos (Bp  -  qp) 


d  LcosA„sinq  LcosA_sinq 

p  _  _ 2  p  _  _ 2  p 

~  sin0p  cos  (Bp  -  <Jp) 


(8.8) 


a  =  —  +  LsinA2tan^c 


(8.9) 


r<»  =  -j-  -  a 
c  2 


LcosA0sinq 
2  r 


-  77  -  LsinA-tan^ 

2  2  c 


(8.10) 


- - — — 7-  +  asini^  ^ 

sin^c  P 


(8.11) 


For  right-hand  gear  we  need  the  following  changes: 


* 

e 

p 


270°  +  Sp  -  qp 


(8.6a) 


b 

P 


LC0SA2 

sin(0* 

P 


cosB 

_ E _ 

-  180°) 


LcosA^ 

cos(8p 


cosB 


-  V 


(8.7a) 


d  LcosA„  sinq 

_E  «  2 _ 

2  sin(0*  -  180® ) 
P 


LcosA^  sinqp 
cos(8p  -  qp) 


(8.8a) 


Let  us  now  consider  point  N  whose  surface  coordinates  of  the  tool  cone 
* 

are  0  and 
P 


(p)  (P)  *  *  * 

u  -  (Th'n  *  (K  'N  +  N  N  =  u  +  Au 
P  P  P 


(8.12) 


Point  N  will  become  the  point  of  contact  if  the  normal  vector  of  this  point 

passes  through  M^.  This  requirement  will  be  observed  if  the  cradle  is  turned 

at  a  certain  angle  4>  that  is  shown  in  Figure  8.3.  We  can  determine  angle  <t> 

P  P 

using  the  equation  of  meshing  of  the  generating  gear  and  the  generated  gear, 
equation  (A. 22),  together  with  equations  (8.5),  (8.6),  (8.7),  (8.8),  (8.9), 
(8.10),  (8.11),  and  (8.12).  This  yields 


f9(Au  ,  <fc  )  * 

2  p  p 

Au  cos(8  +  0  )  ±  LcosA^sin^  ^sin<j>  +  LsinA«cos^  ^  [cosB  -  cos(8  +  4>  )] 
p  p  p  2  Tc  p  2  Tc  p  p  p 

(8.13) 
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IV 


Lsl 


3 


The  upper  sign  corresponds  to  the  left-hand  spiral  bevel  gear  and  the  lower 


to  the  right-hand  spiral  bevel  gear.  Based  on  trigonometry  transformations 


a  closed  form  solution  for  0  is  obtained.  That  is 

P 


^p  -A  i  V  A2  +  B2  -  C2 

tan  *  - 1 - 

2  C  -  R 


Here: 


A  =  cosA„sin'^p)  -  sinA„cosi^  ^sinB  +  sinB 
Z  c  Z  c  p  L  t 


B  =  — r1-  cosB  -  sinA^cosiiJ  ^'cosB 
L  p  2  c  t 


C  *  sinA„cosip  p,cos8 
2  c  F 


For  gear  concave  side  using  the  same  equations  except 


^  -  LsinA„tan(JJ  ^ 
2  2  c 


LcosA0sinq  „  ,  . 

+  a  =*  - — — * - ^  +  -r  -  LsinA9tan4>  p 

siru3*  2  2  c 

P 


we  can  get  the  same  result. 


(8.14) 


(8.15) 


(8.16) 


Equation  (8.14)  determines  the  main  gear  contact  points  for  both  convex 


side  and  concave  side.  In  general,  AUp  >  0  for  gear  convex  side,  and  Au^  <  0 


for  gear  concave  side.  The  absolute  value  of  Au^  for  gear  convex  side  is 


VV 


greater  than  that  of  Au  for  gear  concave  side. 

P 


Direction  of  the  Contact  Path.  The  direction  of  the  contact  path  may  be 
determined  by  the  direction  of  the  tangent  to  the  contact  path  at  the  main  con¬ 
tact  point.  This  can  be  done  by  differentiating  equation  (4.22).  In  the  pro¬ 
cess  of  motion  the  contact  normal  performs  a  parallel  motion.  Thus  parameter 

T  is  a  constant  and  this  yields  that 
P 

dip  =  0  d9p  =  -d<|>  (8.17) 

Taking  into  account  equations  (8.17)  and  (8.5),  we  obtain  after  the  differen¬ 
tiation  of  equation  (4.22)  the  following  equation 

sint  +  b  cos(q  +  4>  )sinii/^  -  b  cost!  ^cosB  tanA_  =  0 
d<$>  p  p  p  p  c  p  c  p  z 


The  upper  sign  corresponds  to  the  left-hand  gear  and  the  lower  to  the  right-hand 

*  * 

P 


gear.  At  the  main  contact  point  we  have  that  0p  =  0p,  cos0  =  sin(8  -  q  ), 


LcosA-cosB 

b  =  _ l _ E_. 

P  cos ( 8  -  q  ) 


Then  we  obtain 


du 

d(J) 


LcosB  r  . 

p  SinTp  +  cos(8p  -  qp)  [cosA2  cos(qp  +  <J>p)sin^P 


sinA0sin(B  -  q  )cosi|/P^l  =  0 
2  P  P  c  J 


(8.18) 


D  =  ~7i~  sinB  sinT„  -k  LcosA„sini{/^cosB  sin<J>  +  LsinA„cos<i/^cos^B 
d<Pp  P  P  2  c  P  P  2  c  p 


UU  /  \  (  X 

E  =  -rr^  cosB  sinx  +  LcosA-simK^'cosB  cosd>  -  LsinA„cos^  ^  sinB  cosB 
d4>p  P  P  2  c  p  p  2  c  P  F 


There  is  a  particular  case  when  the  direction  of  the  contact  path  is 

du  d 

parallel  to  the  root  cone  and  thus  -rr^-  =  0.  For  this  case,  however,  -£ 

d*P  2d 

becomes  substantially  larger  than  LcosA^.  More  favorable  ratio  for  2LcosA" 

du  2 

can  be  obtained  by  decreasing  the  value  of  -rr^-  ,  but  a  more  inclined  path 

p  du 

of  contact  will  occur.  A  reasonable  value  of  -rr-  is  about  -0.5  for  the 

dtp 

P 

left-hand  gear  convex  side,  or  -0.2  for  the  left-hand  gear  concave  side. 


Because  both  sides  of  the  gear  are  cut  by  duplex  method,  the  values  of 
du  du 

-jt^  for  both  sides  are  not  independent,  and  we  have  to  choose  the  value  of  -jr 


for  both  sides  by  a  compromise. 


BASIC  PRINCIPLES  FOR  GENERATION  OF  CONJUGATE  GEAR  TOOTH  SURFACES 


INTRODUCTION 

The  method  for  generation  of  conjugate  spiral  bevel  gear  tooth  surfaces 
is  based  on  the  following  principles: 

(1)  Four  surfaces-two  generating  surfaces  (Zp  and  Zp)  and  pinion  and  gear 
tooth  surfaces  (Z^  and  Z£)  are  In  tangency  at  every  instant.  The  ratio  of 
the  angular  velocities  in  motion  of  surfaces  Zp>  Zp>  Z^  and  Z2  roust  satisfy 
the  following  requirements:  (i)  the  above-mentioned  surfaces  must  be  in 
continuous  tangency  and  (ii)  the  generated  pinion  and  gear  must  transform 
rotation  with  zero  kinematical  errors,  and  (iii)  Point  N  of  contact  of  both 
pairs  of  contacting  surfaces  -  Zp  Zp,  Zj  Z2  “  moves  in  plane  II  (Fig.  9.1) 
in  the  process  of  meshing  and  the  common  normal  to  above-mentioned  surfaces 


performs  parallel  motion  in  plane  II. 

Figure  9.1  shows  the  drawings  that  are  represented  in  the  plane  of 

normals,  II*  Point  A  is  the  point  of  intersection  of  the  gear  head-cutter 

axis  with  plane  II  and  A  is  the  initial  position  of  the  normal  to  surfaces 

Ip  and  I  .  Point  D  is  the  point  of  intersection  of  the  gear  cradle  axis 

with  plane  II.  Simultaneously  point  D  is  the  point  of  intersection  of 

pinion  and  gear  axes.  Axes  of  the  gear  cradle  and  gear  head-cutter  are 

parallel  (the  head-cutter  axis  is  not  tilted  with  respect  to  the  axis  of 

(2) 

the  cradle).  The  cradle  axis,  x  ,  is  perpendicular  to  the  gear  root  cone 

m 

but  it  is  inclined  with  respect  to  the  plane  of  normals,  II. 

Point  C  (Fig.  9.1)  is  the  point  of  intersection  of  the  pinion  tool  cone 
axis  with  plane  II.  Surface  Ip  of  the  pinion  tool  cone  will  be  in  contact 
with  surface  £p  (and  I^)  if  the  common  normal  to  surfaces  Ip  and  I^  passes 
through  point  C.  Simultaneously,  the  pinion  surface,  1^ ,  will  be  in 
contact  with  Ip  (and  Ip  and  at  point  N  if  the  equation  of  meshing  for 
surfaces  1^  and  Ip  is  satisfied  at  N. 

Point  0  that  is  shown  in  Fig.  9.1  is  the  point  of  intersection  of  the 
pinion  cradle  axis  with  plane  II.  Point  0  does  not  coincide  with  point  D  of 
the  gear  cradle  axis.  We  recall  that  specific  machine  tool-settings  for 
the  pinion,  AE^  and  AL^ ,  (Fig.  2.3)  have  to  be  used  to  provide  conjugate 
pinion  and  gear  tooth  surfaces.  This  method  for  generation  provides  the 


collinearity  of  vectors  DO  and  NA.  We  recall  that  the  pinion  cradle  axis, 

X(1\  is  perpendicular  to  the  pinion  root  cone  but  it  is  inclined  with 
m 

(2) 

respect  to  gear  cradle  axis,  ,  and  plane  II.  The  orientation  of  axes 

X^ ^  and  X(2)  with  respect  to  each  other  and  plane  II  depends  on  gear  and 
m  m 

pinion  dedendum  angles,  A^  and  The  pinion  cradle  axis  and  the  pinion 


head-cutter  axis  are  parallel  (the  pinion  cutter  is  not  tilted  with  respect 


to  the  pinion  cradle) 


Figure  9.1  shows  the  instantaneous  positions  of  the  contact  normal  N 
and  points  A  and  C.  In  the  process  of  meshing  the  normal  performs  a 
parallel  motion  in  plane  II  while  point  A  (and  respectively  C)  traces  out  an 
ellipse  with  the  ellipse  symmetry  center  D  (symmetry  center  0, 
respectively).  The  dimensions  and  orientation  of  the  gear  ellipse  centered 
at  D  are  known  since  the  gear  machine  settings  are  given.  The  dimensions 
and  orientation  of  the  ellipse  centered  at  0  must  be  determined  to  provide 
the  desired  parallel  motion  of  the  contact  normal.  Then  it  becomes 
possible  to  determine  the  pinion  machine-tool  settings. 


10.  SATISFACTION  OF  THE  EQUATIONS  OF  MESHING 

Four  surfaces  -  E_,  £  ,  £  and  £_  -  will  be  in  contact  within  the 
t  r  I  2 

neighborhood  of  the  instantaneous  contact  point  if  they  have  a  common 
normal  and  the  following  equations  of  meshing  are  satisfied  at  the  point  of 
contact. 

(I)  The  equation  of  meshing  for  surfaces  and  is  represented  as 
follows  (Fig.  9.1): 


(12)  M  .  (12)  (NK  M  _  r  U2)(n)m1  _  n 
v  •  N  =  (co  xr  )  •  N  =  lu)  r  NJ=U 


(10.1) 


Here:  =  c*/  ^  where  o/ ^  and  a are  the  angular  velocities 

(12) 

of  the  pinion  and  the  gear,  respectively.  Vector  w  represents  the 

(N  ) 

velocity  in  relative  motion.  Vector  r  is  the  position  vector  of  contact 

(12) 

point  N  and  N  is  the  contact  normal.  Vector  ui  must  be  directed  along 
the  pitch  line  to  provide  transformation  of  rotation  with  the  prescribed 
angular  velocity  ratio. 


V - - - 

I 

(2)  The  equation  of  meshing  of  surfaces  Up  and 
equation 


[JP2)r(N>N]  = 


is  represented  by  the 


(10.2) 


(P2)  (P)  (2) 

Equation  (10.2)  is  satisfied  if  ur.  =  oj  -  ur  is  collinear  to  the 
unit  vector  a  of  the  pitch  line.  This  requirement  is  observed  already 
since  the  gear  cutting  ratio  is  determined  with  equation  (8.5). 

(3)  The  equation  of  meshing  of  tool  surfaces  £p  and  1^,  is  represented  as 
follows : 


,/PF) 


N  =  (v<P>  -  v(F>)  •  N  ={(JP)x  r<N))  - 


[(u/F)  X  r(N))  +  (DOxo/F)]  j-  •  N  =  0 


(10.3) 


Deriving  equation  (10.3)  we  substitute  the  sliding  vector  that 

passes  through  0  (Fig.  9.1)  by  an  equal  vector  that  passes  through  D  and 
the  vector  moment  represented  by 


m 


DO  x  u 


(F) 


In  the  process  of  motion  the  contact  normal  keeps  its  original  direction  that 
is  parallel  to  DO.  Thus  equation  (10.3)  yields 


[U<PF>r(B)N]  -  r<JP>  -  „<F>)  rW  m  -  0 


(10.4) 


Equation  (10.4)  may  be  interpreted  kinematically  as  a  requirement  that 
(PF)  (N) 

vectors  oo  ,  r  and  N  must  line  in  the  same  plane,  the  plane  of  normals 


IT.  This  means  that 
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where  i  is  the  unit  vector  of  the  -axis  that  is  perpendicular  to  plan  II. 

It  was  mentioned  above  that  the  proposed  method  for  generation  is  based 
on  the  parallel  motion  of  the  contact  normal  that  slides  along  two  related 
ellipses.  It  was  assumed  that  this  motion  is  performed  with  the  following 
conditions  (see  Section  7): 


p  *  ^  Vi  "  V 


^  .  dp 

dt  dt  dt 


(10.6) 


where  y_  and  u  are  the  ellipse  parameters. 

%  % 

Let  us  prove  that  the  requirement  can  be  observed  if 

|u,(F>|  =  |  and  equation  (10.5)  is  satisfied. 

Vectors  and  are  directed  along  the  -  and  axes  and 

~  mm 

equation  (10.5)  may  be  represented  as  follows 


•  i  .  </p>  • 


(10.7) 


(F)  fP)  (F)  ( 

Since  ur  =  ur  =  oj  we  obtain  that  projections  of  vectors  to  and  ur 

on  the  normal  to  the  plane  are  equal. 

This  yields  that 


^F_^P  dH 
dt  dt  "  dt  ~n 


Here: 


=  to(i(1)  *  i  )  =  o)(i(2)  •  i  ) 
dt  -m  -n  ~m  ~n 


(10.8) 
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Vectors  —  and  represent  the  angular  velocities  of  links  1  and  1*  (Fig. 


7.3)  in  their  rotational  motions.  Equations  (10.8)  also  yields  that 


i  =  i 
~n  ~m 


(10.9) 


Requirement  (10.9)  can  be  satisfied  with  a  specific  orientation  of  the 


pinion  cradle  coordinate  system,  S.  ,  with  respect  to  gear  cradle  coordinate 


system,  Sf.  This  orientation  may  be  achieved  with  a  certain  value  of  angle 
e  (Fig.  5.1,c)  that  may  be  determined  by  using  equation  (10.9).  The  matrix 


representation  of  equation  (10.10)  is  given  as  follows: 


I*  0  °"L„£"Lf,»  1  0 


[1  0  »'  'L„£"Lfh"H» '  '1  0 


(10.10) 


Here,  matrix  [Lf  ]  is  given  by  Eq.  (6.4);  matrix  [L  ]  is  the  transpose 


matrix  of  [L^]  given  by  equation  (6.1);  matrix  [ 1  is  the  sub-matrix  of 


h 1  given  by  equation  (b.l);  matrix  [L^  is  represented  as 


follows  (Fig.  2.3): 


cosAj  0  sinAj 


1  0 


-sinAj  0  cosAj 


(10.11) 


Equation  (10.10)  yields 


AiiH 


(10.12) 


cos  Yj cose  -  sinctanri 


cosA^  -  siny^sin(y^  -  A^) 
cos(Y1  -  A^ 


Equations  (10.12)  provides  two  solutions  for  angle  e  and  the  smaller  value 
of  e  is  to-be  chosen. 


11.  EQUATIONS  OF  SURFACE  TANGENCY  AT  THE  MAIN  CONTACT  POINT 

We  consider  that  the  coordinates  of  the  main  contact  point  N  on  surface 

(p) 

Ip  (Fig.  8.2)  and  the  direction  of  the  surface  unit  normal  n  at  point  N 
are  known.  We  have  to  provide  that  surfaces  Ip,  1^  and  Ip  will  be  in 
tangency  at  the  chosen  point  N.  Surfaces  Ip  and  ^  are  already  in  tangency 
at  N  since  the  equation  of  meshing  of  Ip  and  is  satisifed  at  this  point 
(see  Section  8).  Then  we  have  to  provide  the  tangency  of  surfaces  £p  and 
Ip,  Ip  and  ^  and  ^  and 
Equations  of  Tangency  of  Ip  and  I  p 

The  tangency  of  generating  surfaces  £p  and  Ip  at  point  N  is  provided  if 
the  following  equations  are  satisfied  at  N: 


r(F)(u 


p,  0p,  <J>F,  ^ 


(F) 


.(F) 


*F’ 


AE, 


ALj)  = 


.(P) 


(11.1) 


n(F><v  *F>  *<f>,  V .  „<p> 


(11.2) 


( p)  ( p) 

Here  r  and  n  are  the  given  position  vector  and  surface  Ip  unit 
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normal.  Vector  Eqs.  (11.1)  and  (11.2)  may  be  considered  in  any  coordinate 
system,  for  instance,  in  system  S  Vector  Eq.  (11.1)  provides  three 

independent  scalar  equations  but  Eq .  (11.2)  provides  only  two  ones  since 
Jn |  =  1.  Parameters  Up  and  0p  are  the  surface  Ip  coordinates  for  the  point 
of  tangency.  Parameter  <t>p ,  Is  the  angle  of  cradle  turn  for  the  installment 
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of  the  generating  surface  £p  and  it  may  be  chosen  of  any  value,  Including 

4>p  =*  0.  Parameters  and  bp  determine  the  location  and  orientation  of  the 

pinion  head-cutter  in  coordinate  system  S  (Fig.  2.1);  ij/^  is  the  blade 

ra  c 

angle  (Fig.  2.1);  AE^  and  AL^  are  the  machine-tool  settings  that  have  been 

(F) 

shown  in  Fig.  2.3;  r^  is  the  radius  of  the  head-cutter  circle  obtained  by 

the  intersection  of  the  head-cutter  surface  with  plane  =  0  (Fig.  2.1). 

m 

We  may  represent  vector  equations  of  tangency  (11.1)  and  (11.2)  as 
follows 


■  lK„<1V2)"t<<P2)1  ■  (11.3) 

mm  m 

>”((FlV  -  'L»<1V2)"”<(2))  ■  IL„<1)hllLh£1|LfB(2))["(2)1  <"-4> 

m  m  m 

(2) 

Matrix  [Mc  ]  is  represented  by  the  following  equation  (Fig.  11.1) 
r  m 


cosA^  0  -sinA2  LsinA2COsA2 


-  0 


1  0 


(11.5) 


sinA2  0  cosA2  Lsin 


0  0 


Matrix  [M  ]  is  the  transpose  matrix  of  matrix  [M^l  that  is  represented  by 
equation  (5.1). 

Matrix  [M  ]  is  represented  as  follows  (Fig.  2.3) 
m  n 


The  column  matrix  lr  (2y  and  2  ar6  represented  by  equation  (4.17)  and 
(4.18).  After  matrix  multiplications  we  obtain 

bll  b12  b13  b14 

(2) 1  _  u  k  h  K 

tMm  m  1  b21  b22  b23  b24 

b31  b32  b33  b34 

_0  0  0  1 

Here: 

=  cos£cos(Y1  -  A1)cos(y1  +  A2)  +  sin(y1  -  A1)sin(y1  +  A2) 
b12  =  -sinecos(Yj  -Aj) 

b13  =  -cosccos(y1  "  A1)sin(Y1  +  A2)  +  sin(Y  1  “  A1)cos(y1  +  A2) 

b.  .  =  LslnA.  +  LsinA-  b., 

14  1  2  11 

b21  -  sirecos(Yj  +  A2) 
b22  =  cose 


V  V  •»  •  •'  V  "iV/- V.Y  JV 
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sinesinCYj  +  &>) 


b2^  =  -  AE^+  LsinA2 


b31  =  -cosesin(Y1  -  A1)cos(y1  +  A2>  +  cosCYj  -  A1>sin(Y1  +  A2> 


b32  =  sinesin(Y^  -  A^) 


b  =  cosesin(Y1  “  A^sin^j  +  A2>  +  cos^  -  A^cosCyj  +  A2> 


b_ .  =  -AL,  +  LsinA-  b_, 
34  1  2  31 


Equations  (11.3)  and  (11.4)  are  used  for  the  determination  of  pinion 


machine-tool  settings. 


12.  PINION  MACHINE-TOOL  SETTINGS 

The  pinion  machine-tool  settings  are  represented  by  the  following 
(F) 

parameters:  r^  -  the  radius  of  the  head-cutter  circle  measured  in  plane 

X  =  0;  b„  and  q_  -  parameters  that  determine  the  location  of  the 
m  r  r 

head-cutter  in  plane  X  =0  with  4>  =  0  (see  Fig.  2.1); 


ip  -  parameter  that  determines  the  blade  angle; 


-  the  ratio  by 


cutting;  AE^  and  AL^  -  corrections  of  machine-tool  settings  that  are  shown 
in  Fig.  2.3. 

The  determination  of  the  pinion  machine-tool  settings  is  based  on  the 
equations  that  relate  the  parameters  of  the  gear  and  pinion  ellipses, 
equations  of  tangency  of  the  pinion  and  gear  tooth  surfaces  at  the  main 
contact  point  and  the  equation  of  meshing  by  pinion  cutting.  These 
settings  must  be  determined  for  the  following  4  cases  that  are  represented 
in  the  plane  of  normals,  by  Fig.  12.1  and  Fig.  12.2,  respectively.  Figure 
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> 

V,' 


’0V 

v 

I'ft 

■4 


12.1(a)  corresponds  to  the  case  where  the  gear  is  left-hand,  the  gear  tooth 

contacting  surface  is  convex,  r^  >  r^  and  CN  >  AN.  Figure  12.1(b) 

corresponds  to  the  case  where  the  gear  is  left-hand,  the  gear  tooth 

(  F )  (  p  ) 

contacting  surface  is  concave,  r  <  r  and  CM  <  AN. 

c  c 

Drawings  of  Fig.  12.2  correspond  to  the  cases  where  the  gear  is 

right-hand  and  the  gear  tooth  contacting  surface  is  convex  (Fig.  12. 2, a) 

(  F  )  (  p  ) 

and  concave  (Fig.  12. 2, b).  The  difference  between  r  and  r  with  other 

c  c 

parameters  as  given  determines  the  required  dimensions  of  the  contacting 
ellipse . 

Determination  of  the  Cutting  Ratio  mp 

It  is  mentioned  above  that  the  proposed  method  for  generation  provides: 

(i)  the  simultaneous  tangency  of  four  surfaces:  Ip,  Ip,  1^  and  1^,  where 

Lp  and  Ip  are  the  gear  and  pinion  generating  surfaces  and  Ij  and  are  the 

generated  pinion  and  gear  tooth  surfaces;  (ii)  the  generating  gears 

(  P  )  (  F  ) 

(cradles)  are  rotated  with  the  same  angular  velocities  i.e.  ur  ®  ur 
This  yields  the  following  equation  for  the  pinion  cutting  ratio 


(1)  (2) 

0)  w  m. 


m„-  =  - 7—  (see  equation  8.5) 

P2  cosA^ 


(D  N_ 

w  2 


u/2>  W1 


where  and  Nj  are  the  numbers  of  gear  and  pinion  teeth. 
The  final  expression  for  the  pinion  cutting  ratio  is 


42 


Determination  of  b,, 

- F 

The  procedure  of  computations  of  b  is  based  on  the  following  steps: 

r 

(i)  determination  of  semiaxes  of  the  gear  cradle  ellipse,  ap  and  bp:  (ii) 

determination  of  the  parameter  of  orientation  of  the  gear  cradle  ellipse, 

6p;  (iii)  determination  of  the  parameter  of  orientation  of  the  pinion 

cradle  ellipse,  6  and,  (iv)  determination  of  b  . 

r  F 

Step  1:  Determination  of  semiaxes  of  the  gear  cradle  ellipse. 

The  gear  cradle  ellipse  is  obtained  by  intersection  of  the  cylinder 


of  radius  bp  (Fig.  2.1)  with  the  plane  of  normals.  The  semiminor  axis  of  the 


cradle  ellipse  is  the  cylinder  radius  bp.  The  semimajor  axis  is  represented 


by  the  equation 


3P  = 


i<2)  . 

~m 


(12.2) 


where  the  unit  vector  i 


(2) 


-m 


is  collinear  to  the  cylinder  axis,  i  is  the 

~n 


unit  vector  that  is  perpendicular  to  the  plane  of  normals.  The  scalar 


product  i 


(2) 


m 


•  i^  may  be  represented  as  follows 


i<2)  .  i  = 

'  m  ~n 


[1  0  0][L  (2)  ][Lfn] 
m  r  tn 


(12.3) 


Here,  the  column  matrix  represents  the  unit  vector  i^  in  coordinate  system 


(2)  (2) 

S  and  the  row  matrix  represents  the  unit  vector  i  in  S  . 

n  ~m  ~m 


Equation  (12.3)  yields 
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i 


c^s  cos  n 


(12.4) 


v}‘ 

.'t 


aP  cosA^cosn 


(12.5) 


Step  2:  Determination  of  orientation  of  the  gear  cradle  ellipse. 

o 

We  designate  with  Che  unit  vector  of  the  minor  ellipse  axis  and  with 
5p  the  angle  that  is  formed  by  axis  y^  and  the  major  ellipse  axis  (Fig. 
12.3).  Vectors  b°  and  i^^  are  mutual  perpendicular  and  this  yields  that 


5:. 

I 


1  °  0!ILi2>fHLfnl 


slnSpSinncosA^  +  cos 6  sinA^  =*  0 


(12.6) 


(12.7) 


Step  3:  Determination  of  orientation  of  the  pinion  cradle  ellipse. 


The  unit  vector  of  minor  axis  of  the  pinion  cradle  elllose,  bp,  (Fig.  12.3) 

is  perpendicular  to  the  cvlinder  axis  of  radius  h_.  Thus 

F 


i(1)  •  - 


(12.8) 


1 


The  matrix  representation  of  equation  (12.8)  is  as  follows: 


Wm. 


•  It  ; 


r.v 


V  » 


i 


“  0  Oll'-i"  hI|LhfIILf„>  I  -Sln{F 


(12.9) 


Equation  (12.9)  yields 


(12.10) 


where 


A  =  cosesinyj  -  cosY1tan(y1  ~  A^ 

B  =  cosnsine  +  sinnlcosecosY1  +  sinY^tan(Y1  ~  A j ) ] 


(12.11) 


(12.12) 


Step  4:  Determination  of  bp 

It  is  easy  to  prove  that  the  ratio  between  the  axes  of  both  ellipses  is  the 
same  since 


i  ,  ^  =  i(2)  •  i 

~n  ap  -n  „n 


i(l)  •  i  -  i(2)  •  i  (see  Eq.  (10.9)) 


Thus  : 


cosA„cosn  ’ 


P  cosA,cosn 


(12.13) 


A& 

m 

re* 


The  parameters  of  both  ellipses  have  been  represented  by  equation  (7.19)  as 
follows 


2  .  2 


2  2 


2  .  ,2 


in  (q  -  6p)  +  bpCos  (q  -  6p)  =  d1  + 


Equations  (7.15),  (7.16)  and  (12.12)  yield  that 


22?  ?  sin  (q  6  ) 

dl  +  d2  =  bp  cos  (q  -  5p)  +  - - j- 

cos  A^cos  n 


(12.14) 


(12.15) 


where  q  is  the  angle  that  is  formed  between  vector  DO  and  axis  Yn  (Fig. 
7.2).  We  recall  that  vector  DO  is  collinear  to  the  contact  normal.  Using 
equations  (12.14)  and  (12.15),  we  obtain 


"  2  2  2  2 

cos  (q  -  6p) cos  A^cos  n  +  sin  (q  -  6p) 

3p  2  2  2  2 

cos  (q  -  6  )cos  A„cos  n  +  sin  (q  -  6„) 
r  z  r 


(12.16) 


Determination  of  q_ 
-  F 


The  main  idea  of  determination  of  q_,  is  based  on  identification  in 

r 

plane  ^  =  0  of  projections  of  two  vectors:  OC  (Fig.  12.1  and  Fig.  12.2) 

and  a  that  is  the  unit  vector  of  O^^O^^  (Fig.  12.4).  We  recall  that 
~  ms 

points  0^^  and  0  lie  on  the  axis  of  the  cradle  and  0^^  and  C  lie  on  the 
m  s 

axis  of  the  pinion  head-cutter  (Fig.  12.4).  Points  0  and  C  are  the  points 

of  intersection  of  pinion  cradle  axis,  ^ ,  and  the  pinion  head-cutter 

m 

(F) 

axis,  0g  C,  with  the  plane  of  normals,  IT.  Plane  IT  has  been  shown  in 
Fig.  6.1. 

Unit  vectors  i^\  a  and  OC  lie  in  the  same  plane.  Here  i^^  is  the 
unit  vector  of  axis 


K*;. 

Iv.* 


I 

i 


4  *  I,.1 


[a  i(1)  OCJ  =  a  *  (i(1)  x  OC)  =  0 

-  --m  --  '-m 


(12.17) 


The  determination  of  the  q„  Is  based  on  equation  (12.17)  and  the 

r 

computational  procedure  may  be  presented  as  follows: 

Step  1 :  Representation  of  vector  a  in  coordinate  system 

Figure  12.5(a)  and  Fig.  12.5(b)  show  the  orientation  of  unit  vector  a 
in  plane  in  two  cases  where  the  generating  gear  is  left-hand 

(generated  gear  is  right-hand)  and  generating  gear  right-hand  (generated  is 
left-hand),  respectively.  Vector  a  is  represented  by  the  following  column 
matrix: 


t 


0 


i 

\ 

I 


a 


sinqF 

cosqp 


(12.18) 


Angle  qp  is  measured  clockwise.  This  angle  is  negative  if  measured 
in  opposite  direction  as  shown  in  Fig.  2.1. 

"he  generating  gear  generates  the  member-gear  with  the  same  direction  of  the 
spiral  (the  member-gear  is  down  with  the  respect  to  the  generating  gear  as 
it  is  shown  in  Fig.  2.2).  The  generating  gear  generates  the  pinion  with  the 
opposite  direction  of  the  spir  ll  (the  pinion  is  up  with  the  respect  to  the 
generating  gear  as  it  is  shown  in  Fig.  2.3). 


Step  2:  Representation  of  vector  OC  in  coordinate  system  S  .  Using 

n  ° 

equations  (7,5)  and  (12.13)  we  may  represent  vector  OC  as  follows: 
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I 


OC  -  apCos6pCosypQ  ~  bpSinopSinyp^ 


apSinSpCOS^Q  +  bpcos6FsinyF0 
0 


aF  cos6FcosyF0  -  stn6FsinUFOcosncosA2 


sin6pCosyFQ  +  cos6FsinypQCOsr|cosA2 


(12.19) 


The  sub-script  "FO"  in  ypo  indicates  that  the  initial  position  of  the 
contact  normal  (at  the  main  contact  point)  is  considered. 

Step  3:  Representation  of  vector  OC  in  coordinate  system 

The  coordinate  transformation  in  transition  from  to  ^  is 
represented  by  the  product  of  matrices. 


11 

a12 

a13 

21 

a22 

a23 

31 

a32 

a33 

(12.20) 


Then  vector  OC  may  be  represented  as  follows 


a12b2  +  a13b3 


OC  -  a^2b2  +  a23b3 


(12.21) 


a32b2  +  a33b3 


Step  4:  Representation  of  the  cross-product  im  x  OC. 
It  is  easy  to  verify  that: 


Step  5:  Determination  or  qp 

Equations  '.12.17),  (12.18)  and  (12.22)  yield 


-sinqp(a32b2  +  a33b3)  +  cosqF(a22b2  +  a^b.^)  =  0 


a22b2  +  a23b3 
a32b2  +  a33b3 


(12.23) 


Here : 


a22  =  -sin^cosY^sinn  +  cosCcosH 


(12.24) 


a23  =  -sinesiny^ 


a32  =  sinn [cos£cosy1sin(y1  “  Aj)  -  siny1cos(y1  -  Aj)] 


(12.23) 


r* 

Lv 

Lt»' 

I 


+  cosnsinesin(y  -  A^) 


(12.26) 


a33  =  cosy1cos(y1  -  Aj)  +  cosesiny1sin(y1  -  A^) 


(12.27) 


Expressions  of  b2  and  b3  have  been  represented  by  the  column  matrix  (12.19) 

in  terms  of  Upg*  Our  purpose  is  to  derive  relations  between  Upg»  Dpg  and 

y  and  z  *  Here:  is  the  parameter  of  the  gear  cradle  axis  for 

n  n  rO 


point  A  (Fig.  9.1);  and  z are  coordinates  of  the  main  contact 

point.  Equations  (7.17)  and  (7.18)  yield: 


tan 


—  tan(q  -  6  )  -  cosncosA„ 

d1  r  / 


(12.28 ) 


tan(q  -  6^)  +  —  cosncosA^ 


Here  (see  equations  (7.15)  and  (7.16)) 


d2  tan(q  -  6p)tanupQ  +  cosUcosA2 


dp  tan(q  -  6p)  -  cosr)cosA2tanypQ 


(12.29  ) 


The  determination  of  parameter  ypQ  is  based  on  the  following  considerations 


(Fig.  12.6) 


0  A  =  0  N  +  NA 
n  n 


(12.30) 


Equation  (12.30)  yields 


y^A^  -  y^N^=  -AN  *  cosq 


z(A)  _  z ( N )__  -AN  .  gl 
n  n 


Thus : 


Z(A) 

_z(N) 

n 

n 

(A) 

(N) 

yn 

-yn 

(A)  k7  0 

yn  yn 

Using  equations  (7.4)  and  (12.5)  we  get 


(12.31 ) 


(A)  t  /sinS  c°svp0 

M  c.o3ncosA2  +  cosl5PSlnl,rO 


(12.32) 


(A)  ,  0 Os6PCOSUP0  .  ,  ,  , 

?n  -  bP  \-co?TcosA2  -  S1"6pSlnMPO 


(12.33) 


(N)  (N) 

y  and  z  are  coordinates  of  the  main  contact  point  N  in  coordinate 
n  n  r 

system  S^,q  determines  the  orientation  of  the  contact  normal  in  plane 

x  =  0 . 
n 

Using  Eqs.  (12.31)  -  (12.33)  we  may  determine  parameter  UpQ« 

Determination  of  Relation  Between  AE^  and  AL^ 

Parameters  AE^  and  AL^  of  the  pinion  machine-tool  settings  have  been 

shown  in  Fig.  2.3.  Our  goal  is  to  prove  that  the  proposed  method  for 

generation  relates  AE^  and  AL^  in  a  certain  way. 

Consider  the  drawing  of  Fig.  12.7.  Point  D  coincides  with  the  origins: 

0^  of  coordinate  system  S^,  0^  of  coordinate  system  (Fig.  5.1)  and  the 

origin  of  coordinate  system  S^.  Point  D  is  also  shown  in  Fig.  9.1.  Point 

0  is  the  point  of  intersection  of  the  pinion  cradle  axis  with  the  plane  of 

normals,  II.  Vector  DO  must  be  collinear  to  the  contact  normal  (Fig.  9.1). 

Figure  12.7  shows  a  plane  that  is  drawn  through  points  D  and  0. 

_  m 

Since  vectors  i^\  O^^D  and  0D  lie  in  the  same  plane,  we  have  that 
~m  m 


[ol1 ^D  DO  i  ]  =  0 

m  ~tn 


(12.34) 


Vector  0^  D  is  represented  by  the  column  matrix  (see  Eq.  4.12)) 


LsinA, 


0(  ^D  =  0(  L)0  = 

m  m  h 


(12.35) 


Vector  DO  is  collinear  to  the  unit  contact  normal  and  is  given  by  (see 


Eq.  4.6)) 


r  <F>1 

‘"CD1  ‘ 

m 


(F' 

cos^  ?inT_, 
c  F 


(12.36) 


(F) 

COSljJ  COST- 

c  F 


Equations  (12.34)  -  (12.36)  yield 


LsinA^ 

AEi 

ALj 

sin^(F) 

c 

.(F)  . 

cosip  sini_ 

c  F 

,  (F) 

cosip  CO  ST„ 

c  r 

=  0 

(12.37) 

1 

0 

0 

Equation  (12.36)  results  in  that 


(12.38) 


where 


T  =  0  ±  q  +  <b 

F  F  mF  yF 


(12.39) 


Determination  of  4^  ,  T„  and  AE, 

- rc  j — F - 1 


Parameter  4^  determines  the  blade  angle  of  the  pinion  head-cutter. 


Parameter  T p  with  the  known  value  of  q^  determines  the  cone  surface 


C  i*  7 

parameter  0p.  We  may  determine  4^  and  T  by  using  the  equations  of 
tangency  at  the  main  contact  point  represented  as  follows: 


XT* 


•  *  H"  V5  -  * 


*4‘T' 


II 


V 


v 


| 

$ 


tTl 


v 


ft 


I1 


1 


(12.40) 


(P)  (F) 

r  =  r 


(P)  (F) 

n  =  n 


(12.41) 


Equations  (12.40)  and  (12.41)  provide  the  equality  of  position-vectors 
and  unit  surface  normals,  respectively  at  the  main  contact  point  for  the 
generating  surfaces  Ip  and  Ip.  Vector  equations  (12.40)  and  (12.41) 
provide  three  and  two  scalar  equations,  respectively,  since  In'1  |  =  [  n v  | 
=  l.  Using  vector  Eq.  (12.41)  we  may  determine  parameters  and  ip.  To 

determine  AE^  we  need  only  one  scalar  equation  from  the  three  ones  provided 
by  vector  equations  but  the  remaining  two  scalar  equations  should  be 
checked  to  ensure  that  they  are  satisfied  with  the  determined  machine-tool 
settings . 


13.  INSTALLMENT  OF  MACHINE-TOOL  SETTINGS 

The  installmer t  of  the  eccentric  angle,  EA,  and  cradle  angle,  CA,  on 
the  Gleason  cutting  machine  No.  26  and  No.  116  provides  the  required  values 

of  q.  and  b.  (j  =  P,  F) . 

J  J 

Consider  that  a  left-hand  gear  is  generated.  Figure  13.1  shows  two 
positions  of  the  eccentric  and  its  center:  before  and  after  the 
installment  of  the  cradle  angle,  CA.  The  axis  of  the  cradle  is 
perpendicular  to  the  plane  of  drawings  and  pointed  to  the  observer.  The 
axis  of  the  head-cutter  passes  through  the  center  of  the  cradle  since  the 
eccentric  angle  is  not  installed  yet.  Figure  13.2  shows  two  positions  of 
the  axis  of  the  head-cutter:  before  and  after  the  installment  of  the 
eccentric  angle.  It  is  evident  that  angles  q.,  CA  and  EA  are  related  by  the 
equation 


where  00^  =  k  is  the  given  constant  value.  Figures  13.3  and  13.4 
correspond  to  the  installment  of  a  right-hand  gear  (the  generating  gear  is 
left-hand).  Figure  13.3  shows  the  installment  of  the  cradle  angle  (CA). 
Figure  13.4  shows  the  installment  of  the  eccentric  angle  and  the  relations 
between  angular  parameters  CA,  EA  and  q^ .  Equations  (13.1)  and  (13.2) 
also  for  the  case  of  generation  of  a  right-hand  gear  but  the  value  of  q  is 

j 

larger  than  in  the  case  of  generation  of  a  left-hand  gear  (see  Fig.  13.2). 

The  installment  of  corrections  of  machine-tool  settings  Ae^  and  AL^  are 
only  applied  for  the  generation  of  the  pinion.  These  corrections  are 
installed  on  the  Gleason's  cutting  machine  by  "the  change  of  machine  center 
to  back"  and  "the  sliding  base". 

The  correction  Ae^  represents  the  shortest  distance  between  the  crossed 
axes  of  the  pinion  and  the  gear;  AEj  is  called  "the  offset"  in  Gleason 
terminology  (Fig.  13.5). 

The  correction  Al^  is  installed  at  the  machine  as  a  vector-sum  of:  (i) 
the  change  of  machine  center  to  back  (CB)  and  (ii)  the  sliding  base  (SB) 
(Fig.  I3.6a,b).  The  change  of  machine  center  to  back  is  directed  parallel 
to  the  pinion  axis  and  the  sliding  base  is  directed  parallel  to  the  cradle 
axis.  It  is  evident  that 

AL 

CB  =  - —  ,  SB  =  AL  tany  (13.3) 

cosy  1  R 

K 

Here,  y  is  the  root  cone  angle;  CB  is  the  machine  center  to  back;  SB  is 
K 


the  sliding  base. 


CONCLUSION 


The  authors  proposed  a  method  for  generation  of  spiral  bevel  gears  with 
conjugated  gear  tooth  surfaces  with  the  following  properties: 

(i)  The  transformation  of  rotation  is  performed  with  zero  or  almost 
zero  kinematical  errors. 

(ii)  The  contact  point  of  gear  tooth  surfaces  (the  center  of  instantaneous 
contact  ellipse)  moves  in  a  plane  (it)  of  a  constant  orientation 
that  passes  through  the  pitch  line  of  the  gears. 

(iii)  The  normal  to  the  contacting  surfaces  moves  in  the  process  of  motion 
in  plane  tt  keeping  its  original  orientation. 

(iv)  It  is  expected  that  due  to  the  motion  of  the  contact  ellipse  along 
but  not  across  the  gear  tooth  surfaces  the  contact  ratio  will  be 
increased  and  the  conditions  of  lubrication  improved. 

The  authors  developed  a  method  of  parallel  motion  of  a  straight  line  by 
two  ellipses  with  related  dimensions  and  orientation. 

A  TCA  program  for  the  simulation  of  meshing  and  bearing  contact  for  the 
proposed  gears  has  been  developed.  The  advantage  of  the  proposed  gearing  is 
that  the  gears  can  be  manufactured  by  the  Gleason's  equipment. 
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GEAR  SPIRAL  ON  MACHINE  PLANE 


ELLIPSE 


RIGHT-HAND  GEAR 


IQ)  GEAR  TOOTH  CONVEX  SIDE 


(b)  GEAR  TOOTH  CONCAVE  SIDE 
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HEAD  CUTTER  AXIS 
(  POSITION  I  I 


(  POSITION  2  I 


SPIRAL  DIRECTION 
OF  GENERATING 
GEAR 
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ECCENTRIC 


7m 


NUMERICAL  EXAMPLE 


Basic  Machine  Settings 


16.7979304880 


23.188994098 


Machine  Root  Angle  YR  72.4097056667°  12.0236276667°  12.0236276667° 


3.37751754380  3.31529203349 


Setting  Angle 


73.6837239464°  75.1337128849°  71.9315208492° 


Machine  Offest 


0 . 0050825329 1 A  |  -0 . 0098 1 69  30^ 


Mach.  Ctr.  to  Back 

Xb 

I 

o 

. 

-0.00041258199  !  0.000621-'' 

Sliding  Base 

xs 

0 

-0.00008594705  1  0 .  One  ; 

Ratio  of  Roll 

m 

0.973756061793 

0.237501478486  .  .  2  ‘ 

Orientation  Angle 

.  c 

, 

-0.26906022162' 

Tilt 

i 

0 
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MACHINE  SETTINGS  OF  CONJUGATE  SPIRAL  BEVEL  GEAR 

AUTHORS:  FAYDOR  LITVIN 
WEI-JIUNG  TSUNG 
HONG-TAO  LEE 
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***************************** ********* ******* ********* ******** 


C 

C  PURPOSE 
C 

C  To  find  the  machine  setting  parameters  for  the  generation  of 
C  spiral  bevel  gears. 

C 

C  NOTE 
C 

C  This  program  is  written  in  FORTRAN  77.  It  can  be  compiled  by 
C  V  compiler  in  IBM  mainframe  or  FORTRAN  compiler  in  VAX  system. 
C 

C  DESCRIPTION  OF  INPUT  PARAMETERS 
C 


c. .  . 

JJ 

JJml  for  left-hand  gear,  JJ-2  for  right-hand  gear 

c. .  . 

TNI 

teeth  number  of  pinion 

c. .  . 

TN2 

• 

• 

teeth  number  of  gear 

c. .  . 

D1DG,  D1MIN 

• 

• 

dedendum  angle  of  pinion  (degree  and  arc  minute. 

c. .  . 

respectively) 

c. .  . 

D2DG,  D2MIN 

• 

• 

dedendum  angle  of  gear  (degree  and  arc  minute, 

c. .  . 

respect ively) 

c. .  . 

GAMADG 

• 

• 

shaft  angle  (degree) 

c. .  . 

BPDG 

• 

• 

mean  spiral  angle  (degree) 

c. .  . 

RL 

• 

• 

mean  cone  distance 

c.. . 

W 

• 

• 

point  width  for  gear  finishing 

c. .  . 

ALFAP 

• 

• 

blade  angle  of  gear  cutter  (degree) 

c. .  . 

RCF1 

• 

• 

pinion  cutter  radius  for  pinion  concave  side 

c. .  . 

RCF2 

pinion  cutter  radius  for  pinion  convex  side 

c.. . 

DUP1 

a  chosen  value  for  pinion  concave  side  to  locate  the 

c 

main  contact  point  at  the  desired  location 

c... 

DUP2 

: 

a  chosen  value  for  pinion  convex  side  to  locate  the 

c 

main  contact  point  at  the  desired  location 

c. .  . 

DUPFEP 

: 

a  chosen  value  for  pinion  concave  side  to  get  the 

c 

c 

desired  direction  of  the  contact  path 

C  DESCRIPTION  OF  OUTPUT  PARAMETERS  j 

C  i 


C. .  . 

PHPDG 

:  blade  angle  of  gear  cutter  (degree)  . 

C.  .  . 

PHFDG 

:  blade  angle  of  pinion  cutter  (degree)  j 

C. .  . 

RCP 

:  radius  of  gear  cutter  (measured  on  cradle  plane)  j 

C. .  . 

RCF 

:  radius  of  pinion  cutter  (measured  on  cradle  plane)  r 

C. .  . 

TPDG 

:  the  difference  between  the  sum  of  gear  cone  surface  { 

C 

angle  coordinate  and  generating  surface  rotation 

C 

angle,  and  machine  setting  angle  (degree)  [see 

C 

equation  (4.18)  in  report] 

S 

\ 


nnn  non  non 


C...  TFDG 

C 

C 

C. . . 

C...  BP 
C...  BF 
C...  QPDG 
C...  QFDG 
C...  MP2 
C 

C...  MF1 
C 

C...  DEI 
C...  DL1 
C 

C...  FEEODG 
C 

C...  DELTADG 

C 

C 


:  the  difference  between  the  sum  of  pinion  cone 
surface  angle  coordinate  and  generating  surface 
rotation  angle,  and  machine  setting  angle  (degree) 
[see  equation  (4.3)  in  report] 

:  radial  of  gear  cutter 
:  radial  of  pinion  cutter 
:  machine  setting  angle  of  gear 
:  machine  setting  angle  of  pinion 
:  ratio  of  cradle  angular  velocity  for  cutting  the 
gear  and  gear  angular  velocity 
:  ratio  of  cradle  angular  velocity  for  cutting  the 
pinion  and  pinion  angular  velocity 
:  machine  offset  for  cutting  the  pinion 
:  vector-sum  of  (1)  the  change  of  machine  center  to 
back  and  (2)  the  sliding  base 
:  generating  surface  rotation  angle  at  initial  main 
contact  point 

:  rotation  angle  of  frame  h  relative  to  frame  f  about 
Z1 


IMPLICIT  REAL*8(A-H,0-Z) 


INPUT  THE  DESIGN  DATA 


DATA  JJ/1/ 

DATA  TNI .TN2/10.D00, 41 .D00/ 

DATA  D1DG.D1MIN/1. D00, 41 . D00/ 

DATA  D2DG,D2MIN/3.DOO,53.DOO/ 

DATA  GAMADC , BPDG/90 . D00 , 35 . D00/ 

DATA  RL/3. 226D00/ 

DATA  W/0.08DOO/ 

DATA  ALFAP/20. DOO/ 

DATA  RCF 1/3. 867079296 16D00/ 

DATA  RCF2/4. 05714282715DOO/ 

DATA  DUP1/0. O4278762835870IDO0/ 

DATA  DUP2/-0. 0377637374416016731DOO/ 
DATA  DUPFEP/-0. 526309620098410257DOQ/ 


CONVERT  DEGREES  TO  RADIANS 


CNST-4 . D00*DATAN ( 1 . DOO) / 1 80 . DOO 
Dl- (D1DG+D1MIN/60 . DOO) *CNST 
D2- (D2DG+D2MIN/60.DO0) *CNST 
BP-BPDG*CNST 

CALCULATE  GAMAl  AND  GAMA2 

GAMA-GAMADG*CNST 

RM12-TN2/TN1 

RM21-TN1/TN2 

GMA1-DATAN (DSIN (GAMA) / (RM1 2+DC0S (GAMA) ) ) 
GMA2-DATAN (DSIN (GAMA) / (RM2 1+DC0S (GAMA) ) ) 
IF(GMAl  .LT.  0. ) THEN 
GMA1-180.*CNST+GMA1 
END  IF 


I F (GMA2  .LT.  0 . ) THEN 
GMA2- 180 . *CNST+GMA2 
END  IF 
C 

C  SUBSTITUTE  SIN  AND  COS  FUNCTION  BY  A  SHORT  NAME 
C 

SNBP-DSIN(BP) 

CSBP-DCOS(BP) 

SNDl-DSIN(Dl) 

CSD1-DCOS (Dl) 

SNRlDl-DSIN(GMAl-Dl) 

CSR1D1-DCOS (GMA1-D1) 

SNR1D2-DSIN (GMA1+D2) 

CSR1D2-DCOS (GMA1+D2) 

SNR2D2-DSIN (GMA2-D2) 

CSR2D2-DCOS (GMA2-D2) 

CSD2-DCOS(D2) 

SND2-DSIN(D2) 

CSGM1-DCOS (GMA1) 

SNGMl-DSIN(GMAl) 

CSGM2«DCOS(GMA2) 

SNGM2«DSIN(GMA2) 

C 

C  CALCULATE  CUTTING  RATIOS 
C 

RMP2-SNGM2/CSD2 
RMFl-TNl*SNGM2/ (TN2*CSD2) 

C 

C  1 1- 1 :  THE  PINION  CONCAVE  PART  ANALYSIS,  11*2:  THE  PINION  CONVEX  PART 
C  ANALYSIS 
C 

DO  10000  11-1,2 
IF (I I  .EQ.  1)THEN 
RCF-RCF1 
DUP-DUP1 
ELSE 

RCF-RCF2 
DUP-DUP2 
END  IF 
C 

C  CALCULATE  PHP 
C 

I F ( 1 1  .EQ  .1)THEN 
PHP-ALFAP*CNST 
ELSE 

PHP- ( 180 . DOO-ALFAP) *CNST 
END  IF 

SNPHP-DSIN (PHP) 

CSPHP-DCOS (PHP) 

C 

C  CALCULATE  FEEP  AT  INITIAL  MAIN  CONTACT  POINT 
C 

A-CSD2*SNPHP-SND2*CSPHP*SNBP+DUP*SNBP/RL 
IF ( J J  .EQ.  2) THEN 
A— A 
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ooo  nno  noo 


END  IF 

B-DUP*CSBP/RL-SND2*CSPHP*CSBP 

C«SND2*CSPHP*CSBP 

D-DSQRT (A*A+B*B-C*C) 

E-C-B 

F1«DATAN((-A+D)/E) 

F2-DATAN((-A-D)/E) 

IF (DABS (FI)  .LT.  DABS (F2) ) THEN 
FEE-2. DOO*F1 
ELSE 

FEE-2. D00*F2 
END  IF 

CSFEE-DCOS (FEE) 

SNFEE-DSIN (FEE) 

CALCULATE  TP 

IF (J J  .EQ.  1)THEN 
TP-90. D00*CNST-BP+FEE 
ELSE 

TP-270. D00*CNST+BP+FEE 
END  IF 

SNTP-DSIN (TP) 

CSTP-DCOS (TP) 

CALCULATE  QP 

IF (II  .EQ.  1)THEN 
IF ( J J  .EQ.  1 ) THEN 

D-DUPFEP*SNBP*SNTP+RL*CSD2*SNPHP*CSBP*SNFEE+RL*SND2*CSPHP 
#  *CSBP*CSBP 

ELSE 

D-DUPFEP*SNBP*SNTP-RL*CSD2*SNPHP*CSBP*SNFEE+RL*SND2*CSPHP 
//  *CSBP*CSBP 

END  IF 

E-DUPFEP*CSBP*SNTP+RL*CSD2*SNPHP*CSBP*CSFEE-RL*SND2*CSPHP 
if  *SNBP*CSBP 
QP-DATAN (-E/D) 

IF(QP  .LT.  0.D00) THEN 
QP- 180 . D00*CNST+QP 
END  IF 
END  IF 

IF(J J  .EQ.  1)THEN 
SNQPFE-DSIN (QP-FEE) 

CSQPFE-DCOS (QP-FEE) 

THP-90.D00*CNST+QP-BP 

ELSE 

SNQPFE-DSIN (QP+FEE) 

CSQPFE-DCOS (QP+FEE) 

THP-270. D00*CNST-QP+BP 
END  IF 

SNTHP-DSIN (THP) 

CALCULATE  RCP 
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CM0«RL*CSD2*DSIN (QP) /DCOS (BP-QP) 

IF (I I  .EQ.  1)THEN 

RCP=CM0-W/2 . 0D00-RL*SND2*SNPHP/CSPHP 
ELSE 

RCP=CM0+W/2 . D00-RL*SND2*SNPHP/CSPHP 
END  IF 
C 

C  CALCULATE  RBP 
C 

RBP=RL*CSD2*CSBP/DCOS (BP-QP) 

C 

C  CALCULATE  UP 
C 

UP=  (RCP*CSPHP/SNPHP+RL*SND2) *CSPHP+RBP* (SNPHP*SNQPFE-CSPHP* 
#  SNTHP*SND2/CSD2) /SNTP 
C 

C  CALCULATE  RN 
C 

AN— SNPHP*CSD2+CSPHP*CSTP*SND2 

BN«CSPHP*SNTP 

RN=DATAN (AN/BN) 

SNRN-DSIN(RN) 

CSRN-DCOS(RN) 

C 

C  CALCULATE  DLTA 
C 

A-CSGM1 

B—SNRN/CSRN 

C“ (SNGM1*SNR1D1-CSD2) /CSR1D1 
XI—  B+DSQRT  (A**2+B**2~C**2) 

X2— B-DSQRT  (A**2+B**2-C**2) 

Y-C-A 

AE"=X1/Y 

BE-X2/Y 

DLT2A-DATAN (AE) 

DLT2B-DATAN (BE) 

IF (DABS (DLT2A) . LT . DABS (DLT2B) ) THEN 
DLTA-DLT2A*2.D00 
ELSE 

DLTA-DLT2B*2.DOO 
END  IF 

SNDLTA-DSIN(DLTA) 

CSDLTA-DCOS (DLTA) 

C 

C  CALCULATE  DETAP 
C 

A01-SND2/CSD2 
B01— SNRN 

DETAP-DATAN (AO 1 /BO 1 ) 

C 

C  CALCULATE  DTAF 
C 


A2I1-CSGMI*SNR1D1/CSR1D1-CSDLTA*SNGM1 

B211— (SNDLTA*CSRN+SNRN*(CSDLTA*CSGM1  +  SNGM1’vSNR1D1/CSR1D1)) 
DETAF-DATAN (a21 I/B21 1) 


non  non  onnonn  nno  n  non 


CSDTAF=DCOS (DETAF) 

SNDTAF=DSIN (DETAF) 

C 

C  CALCULATE  Q 
C 

C  calculate  normal  vector  at  main  contact  point  in  F  coordinate  system 
C 

RNXF2=SNPHP*CSD2-CSPHP*CSTP*SND2 
RNYF2=CSPHP*SNTP 

RNZF2-SNPHP*SND2+CSPHP*CSTP*CSD2 

calculate  normal  vector  at  main  contact  point  in  N  coordinate  system 

RNYN2—SNRN*RNXF2+CSRN*RNYF2 
RNZN2=RNZF2 

Q=DATAN2 (-RNZN2 , -RNYN2) 

QAP=Q-DETAP 
QAF=Q-DETAF 
SNQAP=DSIN (QAP) 

CSQAP=DCOS (QAP) 

SNQAF=DSIN (QAF) 

CSQAF=DCOS (QAF) 

CALCULATE  BF 

BF=DSQRT ( (CSQAP**2*CSD2**2*CSRN**2+SNQAP**2) 

//  /  (CSQAF**2*CSD2**2*CSRN**2+SNQAF**2) )  *RBP 

CALCULATE  PO 

RAP=RBP/(CSD2*CSRN) 


calculute  coordinate  of  main  contact  point  in  M  coordinate  system 

RMXM2=RCP*CSPHP/SNPHP-UP*CSPHP 

RMYM2»UP*SNPHP*SNTP-RBP*SNQPFE 

RMZM2»UP*SNPHP*CSTP+RBP*CSQPFE 

calculute  coordinate  of  main  contact  point  in  F  coordinate  system 

RMXF2-CSD2*RMXM2-SND2*RMZM2+RL*SND2*CSD2 

RMYF2-RMYM2 

RMZF2«SND2*RMXM2+CSD2*RMZM2+RL*SND2*SND2 


calculute  coordinate  of  main  contact  point  in  N  coordinate  system 

RMXN2-CSRN*RMXF2+SNRN*RMYF2 
RMYN2— SNRN*RMXF2+CSRN*RMYF2 
RMZN2-RMZF2 

AG— RAP*SNQAP 
BG-RBP*CSQAP 

CG-RMYN2*DSIN (Q) -RMZN2*DC0S (Q) 

G1 1— BG+DSQRT (AG**2+BG**2~CG**2) 


G22“-BG-DSQRT(AG**2+BG**2-CG**2) 

G33=CG-AG 
P01-DATAN (G1 1/G33) 

P02-DATAN (G22/G33) 

I F (DABS (P01) .LT. DABS (P02) ) THEN 
P0“P02*2.D00 
ELSE 

P0=P01*2.D00 
END  IF 

SNPO=DSIN(PO) 

CSPO-DCOS(PO) 

C 

C  CALCULATE  FO 
C 

AF=BF/ (CSD2*CSRN) 

RNF=BF**2 

RMF-AF**2 

DDll=RAP*SNQAP*CSPO-RBP*CSQAP*SNPO 
DD22=RAP*SNQAP*SNP0+RBP*CSQAP*CSP0 
SSS= (AF*SNQAF) **2+ (BF*CSQAF) **2 
CSFO= (AF*DD11*SNQAF+BF*DD22*CSQAF)/SSS 
SNFO-  (AF*DD22*SNQAF-BF*DD1  1*CSQAF)  /SSS 
F0=DATAN2 (SNFO, CSFO) 

C 

C  CALCULATE  QF 
C 

AA22— SNDLTA*CSGMl*SNRN+CSDLTA*CSRN 
AA23=-SNDLTA*SNGMl 

AA32*SNRN* (CSDLTA*CSGM1 *SNR 1D1-SNGM1 *CSR 1 D1 ) +CSRN*SNDLTA*SNR 1D1 
AA33=CSGM1*CSR1D1+CSDLTA*SNGM1*SNR1D1 
BB22-AF* (CSDTAF*CSF0“SNDTAF*SNF0*CSRN*CSD2) 
BB33-AF*(SNDTAF*CSF0+CSDTAF*SNF0*CSRN*CSD2) 

UUU—  (AA22*BB22+AA23*BB33) 

DDD-AA32*BB22+AA33*BB33 
QFFE-DATAN2 (UUU , DDD) 

QF-QFFE 

SNQF-DSIN(QF) 

CSQF-DCOS(QF) 

SNQFFE-DSIN(QFFE) 

CSQFFE-DCOS (QFFE) 

C 

C  calculate  normal  vector  at  main  contact  point  in  H  coordinate  system 
C 

RNXH2-(CSDLTA*CSGM1*CSGM1+SNGM1*SNGM1)*RNXF2-SNDLTA*CSGM1*RNYF2 
#  +CSGM1*SNGM1*(1 ,D00-CSDLTA)*RNZF2 

RNYH2-SNDLTA*CSGMl*RNXF2+CSDLTA*RNYF2-SNDLTA*SNGMl*RNZF2 
RNZH2-CSGM1*SNGM1*(1.DOO-CSDLTA)*RNXF2+SNDLTA*SNGM1*RNYF2+ 
if  (CSDLTA*SNGMl*SNGMl+CSGMl*CSGMl)*RNZF2 

C 

C  CALCULATE  PHF 
C 

PHF«DASIN(CSD1*RNXH2-SND1*RNZH2) 

IF (I  I  .EQ.  2) THEN 
PHF- 180 . DOO*CNST-PHF 
END  IF 


nnn  onon  onnon  nnn  nnn 


SNPHF=DSIN (PHF) 
CSPHF=DCOS (PHF) 


CALCULATE  TF 

IF(II  .EQ.  2) THEN 

TF=DATAN2 (-RNYH2 , -SND1 *RNXH2-CSDl *RNZH2) 
ELSE 

TF=DATAN2 (RN YH2 , SND1 *RNXH2+CSD1 *RNZH2) 
END  IF 

SNTF=DSIN(TF) 

CSTF=DCOS(TF) 

CALCULATE  THF 

THF=TF+QFFE 

SNTHF=DSIN(THF) 

CSTHF=DCOS (THF) 

CALCULATE  UF 


calculute  coordinate  of  main  contact  point  in  H  coordinate  system 

RMXH2=(CSDLTA*CSGM1*CSGM1+SNGM1*SNGM1)*RMXF2-SNDLTA*CSGMPVRMYF2 
//  +CSGM1*SNGM1*  (1 .  DOO-CSDLTA)  *RMZF2 

RMYH2=SNDLTA*CSGM1*RMXF2+CSDLTA*RMYF2-SNDLTA*SNGM1*RMZF2 
RMZH2=CSGM1*SNGM1*(1. DOO-CSDLTA)  VtRMXF2+SNDLTA*SNGMl*RMYF2+ 

//  (CSDLTA*SNGMl*SNGMl+CSGMl*CSGMl) *RMZF2 

UF= (RCF*CSPHF/SNPHF-RMXH2*CSD1+RMZH2*SND1-RL*SND1) /CSPHF 

CALCULATE  DE1.DL1 

DEl=RMYH2-UF*SNPHF*SNTF+BF*SNQFFE 

DL1=DE1*CSTF/SNTF 

CONVERT  RADIANS  TO  DEGREES 

PHPDG'PHP/CNST 

TPDG=TP/CNST 

QPDG=QP/CNST 

RNDG=RN/CNST 

DLTADG*=DLTA/ CNST 

DTAPDG=DETAP/CNST 

DTAFDG=DETAF/CNST 

QDG'Q/CNST 

PODG'PO/CNST 

F0DG=F0/CNST 

QFDG'QF/CNST 

PHFDG-PHF/CNST 

TFDG-TF/CNST 

FEEDG-FEE/CNST 

WRITE (6, 60000) 

WRITE(6, 60001) 

WRITE(6, 60002) 


102 


IF(II.EQ.I)  THEN 
WRITE(6, 60003) 

ELSE 

WRITE (6, 60004) 

END  IF 

WRITE (6,60002) 

WRITE(6, 60001) 

WRITE (6,60000) 

WRITE(6, 90001) 

WRITE (6 , 90002) PHPDG , PHFDG , RCP , RCF , TPDG , TFDG 

WR I TE (6 , 90003) RBP , BF, QPDC , QFDG, RMP2 , RMF1 , DE 1 , DL 1 , FEEDG , DLTADG 
WRITE (6 , 1) RNDG 
WRITE(6, 20001) 

WRITE (6 , 20002) DTAPDG , PODG , RBP , RAP 
WRITE (6, 20003) 

10000  WRITE (6 , 20004) DTAFDG , FODG , BF , AF 

90003  FORMAT (1H  , 'BP  = ' ,G18 . 12 , 15X, ' BF  = ' , G18 . 12 , / 

#  1H  ,'Q PDG  =' ,G18. 12, 15X, ’QFDG  =',G18.12,/ 

//  1H/MP2  ='  ,G18. 12, 15X,  'MF1  -'.G18.12,/ 

//  1H  ,  '  DEI  =  '  ,  G18 . 12 , 15X, '  DLl  =’,G18.12,/ 

#  1H  , 'FEEODG3' ,G18. 12.15X, 'DELTA  =' ,G18. 12/) 

1  FORMAT (1H  ,' ROTATION  ANGLE  OF  RN  =',G18.12/) 

20001  FORMAT (1H  ,'  GEAR  CRADLE  ELLIPSE  PARAMETER  ' /) 

20002  FORMAT (1H  , ' GEAR  ELLIPSE  ORIENTATION3 ' , G24 . 17 , / 

#  1H  , ' INITIAL  POSITION3 ' , G24 .11,/ 

//  1H  , 'MINOR  AXIS3' ,G24. 17, 5X, 'MAJOR  AXIS3  ’  ,  G24 . 17/) 

20003  FORMAT (1H  ,'  PINION  CRADLE  ELLIPSE  PARAMETER  ' /) 

20004  FORMAT (1H  , 'PINION  ELLIPSE  ORIENTATION= ' ,G24. 17, / 

//  1H, 'INITIAL  POSITION3' ,G24. 17,/ 

//  1H  ,  'MI NOR  AXIS3’ ,G24. 17, 5X, 'MAJOR  AXIS3  ’  , G24 . 1 7/) 

60000  FORMAT (1H1 

60001  FORMAT (1H 

60002  FORMAT (1H 

60003  FORMAT (1H 

60004  FORMAT (1H 

90001  FORMAT (1H 

90002  FORMAT (1H 

//  1H 

#  1H 
END 


»  it  *  *  it  it  *  *  it  it  it  *  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  it  »  ) 
»  it  it  «  ) 

'*  CONCAVE  PINION  PART  ANALYSIS  *') 

'*  CONVEX  PINION  PART  ANALYSIS  *') 

'GEAR  PARAMETER' ,27X, 'PINION  PARAMETER’) 

'PHPDG  =' ,G18. 12, 15X, 'PHFDG  3',G18.12,/ 

'RCP  3 ' , G18 . 12 , 15X, ' RCF  =',G18.12,/ 

'TPDG  =' ,G18. 12, 15X, 'TFDG  =',G18.12) 
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C 

C  PURPOSE 
C 

C  Tooth  contact  analysis  of  spiral  bevel  gears 

C 

C  NOTE 
C 

C  This  program  is  written  in  FORTRAN  77.  It  can  be  compiled  by 
C  V  compiler  in  IBM  mainframe. 

C 

C  This  program  calls  ZSPOW,  a  subroutine  of  IMSL  package. 

C 

C  DESCRIPTION  OF  INPUT  PARAMETERS 
C 


c. . . 

JJ 

: 

JJ=1  for  left-hand  gear,  JJ=2  for  right-hand  gear 

c. . . 

II 

: 

11=1  for  pinion  concave  side,  11=2  for  pinion  convex 

c 

side 

c. .  . 

TNI 

teeth  number  of  pinion 

c. . . 

TN2 

• 

teeth  number  of  gear 

c. .  . 

DlDG,  DlMIN 

• 

• 

dedendum  angle  of  pinion  (degree  and  arc  minute. 

c. . . 

respectively) 

c. . . 

D2DG,  D2MIN 

: 

dedendum  angle  of  gear  (degree  and  arc  minute. 

c. . . 

respectively) 

c. . . 

GAMADG 

shaft  angle  (degree) 

c. . . 

RL 

mean  cone  distance 

c. . . 

PHPDG 

blade  angle  of  gear  cutter  (degree) 

c. . . 

PHFDG 

blade  angle  of  pinion  cutter  (degree) 

c. . . 

RCP 

radius  of  gear  cutter  (measured  on  cradle  plane) 

c. . . 

RCF 

radius  of  pinion  cutter  (measured  on  cradle  plane) 

c. .  . 

TPDG 

the  difference  between  the  sum  of  gear  cone  surface 

c 

angle  coordinate  and  generating  surface  rotation 

c 

angle,  and  machine  setting  angle  (degree)  [see 

c 

equation  (A. 18)  in  report] 

c. . . 

TFDG 

• 

the  difference  between  the  sum  of  pinion  cone 

c 

surface  angle  coordinate  and  generating  surface 

c 

rotation  angle,  and  machine  setting  angle  (degree) 

c. .  . 

[see  equation  (A. 3)  in  report] 

c. . . 

BP 

radial  of  gear  cutter 

c. . . 

BF 

radial  of  pinion  cutter 

c. . . 

QPDG 

machine  setting  angle  of  gear 

c. .  . 

QFDG 

• 

machine  setting  angle  of  pinion  [ALWAYS  consider  it 

c 

as  a  POSITIVE  input  even  it  is  negative  gotten  from 

c 

the  output  of  maching  setting  program 

c. . . 

MP2 

« 

• 

ratio  of  cradle  angular  velocity  for  cutting  the 

c 

gear  and  gear  angular  velocity 

I 


non 


c. .  . 
c 

MF1 

: 

ratio  of  cradle  angular  velocity  for  cutting  the 
pinion  and  pinion  angular  velocity 

c. .  . 

DEI 

: 

machine  offset  for  cutting  the  pinion 

c. . . 
c 

DL1 

• 

vector-sum  of  (1)  the  change  of  machine  center  to 
back  and  (2)  the  sliding  base 

c. . . 
c 

FEEODG 

: 

generating  surface  rotation  angle  at  initial  main 
contact  point 

c. . . 
c 

DELTADG 

rotation  angle  of  frame  h  relative  to  frame  f  about 
Z1 

c. .  . 

DFEPDG 

: 

increment  of  FEPDG 

C 

IMPLICIT  REAL*8(A-H,0-Z) 

INTEGER  NSIG, IER , ITMAX, M, MLOOP , N1 , N2 , N5 , NSIG1 
REAL *8  PAR5 (5) , X (5) , FN0RM5 , WK5 (75) ,MP2,MFl 
EXTERNAL  FCN1 , FCN2 
COMMON/Al/CNST 

COMMON/ A2 /DEI , DL1 , RCP , RL , RCF, PHFDG 

COMMON/A3/SNPHF , CSPHF , SNR 1D1 , CSR 1D1 , SNBP , CSBP , SNR 1D2 , CSR 1D2 
C0MM0N/A4/W , ECENDG , DISTEC 

COMMON /A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/ A12/UF , BF , RBP , UP , RMFl , RMP2 
COMMON/ A14/SNRN , CSRN , RNDG.RN 

C0MM0N/A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A27/XF1 , YF1 ,ZF1 ,XF2, YF2.ZF2 

COMMON/A32/AYN , AZN , BYN , BZN , CYN , CZN , DYN , DZN , RMYN , RMZN 
COMMON/A33/PHP1 , PHF1 , PHP2 , PHF2 , PHP1DG, PHF1DG, PHP2DG , PHF2DG 
COMMON / A3  4 / FE 1 , SNFE1 , CSFE 1 , FE2 , SNFE2 , CSFE2 , FEP , FE22 1 , FE 1 1 1 , 

#  FEPDG,FEFDG,F1UDG,F221DG 

COMMON /A3 5 /SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
COMMON/ A36/RK 1 2 , RK22 , RK1 1 , RK2 1 , UTX 1 1 , UTX 1 2 , UT Y 1 1 , UT Y 1 2 , UTZ 1 1 , 

//  UTZ  1 2 ,  UTX2 1 ,  UTX22 ,  UT  Y2 1 ,  UT  Y22 ,  UTZ2 1 ,  UTZ22 

COMMON/A37/RNXF2 , RNYF2 , RNZF2 , RNXF1 , RNYF1 , RNZF1 
COMMON/A38/DEF, SIGMDG , ALPHDG , AX ISA, AX I SB , ERROR 
COMMON/A40/GMENRA, PMENR1 , PMENR2 
COMMON / A49 /SNQAP .CSQAP , SNQAF ,CSQAF 

COMMON/A55 /DETAP , DETAF , SNDTAP , CSDTAP , SNDTAF , CSDTAF , DTAPDG , DTAFDG 
COMMON / A7  7 /Q , QDG , PODG , FODG , PMIN , PMAX , FMIN, FMAX 

INPUT  DATA 

JJ-1 
II-l 

TNl-lO.DOO 
TN2-41.DOO 
D1DG-1.DOO 
D1MIN-41.D00 
D2DG-3.DOO 
D2MIN-53.DOO 
GAMADG-90 . DOO 
RL-3.226D00 
PHPDG-20.0D00 
PHFDG- 16 . 7979304880D0O 
RCP-3 .83760775903DO0 


c 

C  END 
C 


RCF-3.86707929616D00 
TPDG=53. 0019620800D00 
TFDG=5 1.55 1973 1416D00 
BP=3. 37751 754380DOO 
BF=3. 31529203349D00 
QPDG-73 . 6837239464DOO 
QFDG=77. 1317508049D00 
MP2-.973756061793D00 
MFl= . 237501478486D00 
DE1=0. 508253291408D-02 
DL1=0.403530720434D-02 
FEP0DG=-1 .9980379199559D00 
DLTADG=- . 269060223623D00 
DFEPDG* 1 . OODOO 

OF  INPUT  DATA 

DEF*=0 . 00025D00 
N5-5 

NSIG1«12 

ITMAX=200 

RBP=BP 

RMP2=MP2 

RMF1=MF1 

CONVERT  DEGREE  TO  RADIANS 


CNST=4 . DOO*DATAN ( 1 . DOO) / 1 80 . DOO 

GAMA“GAMADG*CNST 

QP”QPDG*CNST 

PHP-PHPDG*CNST 

PHF=PHFDG*CNST 

TP=TPDG*CNST 

TF=TFDG*CNST 

DLTA“DLTADG*CNST 

QF-QFDG*CNST 

FEPO“FEPODG*CNST 

FEFO-O.DOO 

Dl= (DlDG+DlMIN/60 . DOO) *CNST 
D2= (D2DG+D2MIN/60. DOO) *CNST 
RM12-TN2/TN1 
RM21-TN1/TN2 

GMAl-DATAN (DSIN (GAMA) / (RM12+DC0S (GAMA) ) ) 
GMA2'DATA N (DSIN (GAMA) / (RM2 1 +DCOS (GAMA) ) ) 

M- (360 . DOO/ (2 . 0D00*TN2) ) /DFEPDG 
CALL  HEAD1(JJ,II) 

SUBSTITUTE  SIN  AND  COS  FUNCTION  BY  A  SHORT  NAME 

SNDl-DSIN(Dl) 

CSDl-DCOS(Dl) 

SNRlDl-DSIN(GMAl-Dl) 

CSR1D1-DC0S (GMA1-D1) 

SNR1D2“DSIN (GMA1+D2) 

CSR 1 D2-DC0S (GMA1+D2) 


n  n  n  n 


SNR2D2-DSIN (GMA2-D2) 

CSR2D2-DCOS (GMA2-D2) 

CSD2«DCOS(D2) 

SND2»DSIN(D2) 

CSGM1-DC0S (GMAl) 

SNGMl-DSIN(GMAl) 

CSGM2-DCOS(GMA2) 

SNGM2-DSIN(GMA2) 

SNPHP=DSIN (PHP) 

CSPHP=DCOS(PHP) 

SNQP-DSIN(QP) 

CSQP-DCOS (QP) 

SNTP=DS IN (TP) 

CSTP=DCOS (TP) 

SNDLTA=DSIN (DLTA) 

CSDLTA-DCOS (DLTA) 

SNQF-DSIN(QF) 

CSQF-DCOS(QF) 

SNPHF=DSIN (PHF) 

CSPHF“DCOS (PHF) 

SNTF-DSIN(TF) 

CSTF“DCOS (TF) 

TOOTH  CONTACT  ANALYSIS 

DETERMINATION  OF  THE  KINEMATICAL  ERROR,  BEARING  CONTACT 

TP0=TP 
TFO-TF 
QPO-QP 
QFO=QF 

RINITP*M*DFEPDG 
FEP=-RINITP*CNST+FEPO 
C...  INITIAL  GUASS  FOR  ZSPOW 
X(l)-TFO 
X  (2) *TPO 

X (3) — RINITP*CNST+FEFO 

X(A)«O.ODOO*CNST 

X (5) “0 . 0D00*CNST 

MLOOP*=M*2+ 1 

MLOOPO-M+1 

DO  99  J-l.MLOOPC 

IF ( JJ  .EQ.  1)THEN 

CALL  ZSPOW (FCN 1 , NS I G 1 , N5 , I TMAX , P AR5 , X , FNORM5 , WK5 , IER) 
ELSE 

CALL  ZSPOW (FCN2,NSIG1 , N5 , I TMAX, PAR5 , X, FNORM5 , WK5 , IER) 
END  IF 

CALL  ANGLE (X ( 1 ) ) 

CALL  ANGLE (X (2)) 

CALL  ANGLE (X (3)) 

CALL  ANGLE (X (A)) 

CALL  ANGLE (X (5)) 

FEF-X(3) 

FE2* (FEP-FEPO) /RMP2 
FEl* (FEF-FEFO) /RMFI 
FElll-X(A) 


FE221-XC5) 

FE2PR-FE2-FE221 

FE1PR=FE1~FE1 1 1 

FEPDG-FEP/CNST 

FE1PRD-FE1PR/CNST 

FE2PRD-FE2PR/CNST 

FEP-FEP+DFEPDG*CNST 

CONTINUE 

FE1PR0-FE1PR 

FE2PR0-FE2PR 

RINITP=M*DFEPDG 

FEP=-RINITP*CNST+FEPO 

QPO-QPO 

QFO-QFO 

INITIAL  GUASS  FOR  ZSPOW 
X (1) “TFO 
X (2) -TPO 

X(3)—RINITP*CNST+FEFO 
X  (4)  “0 .  ODOO*CNST 
X(5)-0.0D00*CNST 
MLOOP-M*2+ 1 
DO  88  J* 1 , MLOOP 
IF(J J  .EQ.  1)THEN 

CALL  ZSPOW  (  FCN 1 ,  NS  I G1 ,  N5 ,  I  TMAX ,  PAR5  ,  X ,  FNORM5  ,  WK5  ,  IER) 

ELSE 

CALL  ZSPOW (FCN2 , NSIG1 , N5 , ITMAX, PAR5 , X , FNORM5 , WK5 , I ER) 

END  IF 

CALL  ANGLE (X ( 1) ) 

CALL  ANGLE (X (2)) 

CALL  ANGLE (X (3)) 

CALL  ANGLE (X  (A) ) 

CALL  ANGLE (X (5)) 

FEF-X (3) 

FE2= (FEP-FEPO) /RNP2 
FE1'(FEF-FEF0)/RMF1 
SNFEl'DSIN(FEl) 

CSFEl-DCOS (FE1) 

SNFE2=DSIN (FE2) 

CSFE2-DCOS(FE2) 

SNTF-DSIN (X ( 1) ) 

CSTF-DCOS (X(l)) 

SNTP-DSIN  (X  (2)  ) 

CSTP”DCOS (X (2) ) 

FElll-X(A) 

FE221-X (5) 

FE2PR-FE2-FE221 

FE1PR=FE1-FE111 

ERROR" (FE2PR *3600. DOO-FE2PR 0*3600. D00- (FE1PR*3600. D00-FE1PR0* 
#  3600. D00) *TN1/TN2) /CNST 

TFDG-X(1)/CNST 
TPDG*X (2) /CNST 
FEFDG”X(3) /CNST 
FEPDG-FEP/CNST 
FI  1 1DG-FE1 1 1/CNST 
F221DG-FE221/CNST 


non  non  non 


FE1PRD-FE1PR/CNST 

FE2PRD-FE2PR/CNST 

,  .  DETERMINATION  OF  THE  BEARING  CONTACT 

CALL  FNORM (DN 1 XTH , DN 1XFE , DN 1 YTH , DN 1 YFE , DN 1ZTH , DN 1 ZFE) 

CALL  PNORM (DN2XTH , DN2XFE , DN2 YTH , DN2 YFE , DN2ZTH , DN2ZFE) 

CALL  DUP (UPP , DUPTHP , DUPFEP) 

CALL  DUF (UFF,DUFTHF,DUFFEF) 

CALL  DR 2 (UPP , DUPTHP , DUPFEP , DX2THP , DX2FEP , DY2THP , DY2FEP , 

//  DZ2THP ,  DZ2FEP) 

CALL  DR  1 (UFF , DUFTHF , DUFFEF , DX1THF , DX1 FEF , DY 1 THF , DY 1 FEF , 

//  DZ1THF.DZ1FEF) 

CALL  PC (DN2XTH , DN2XFE , DX2THP , DX2FEP , DN2YTH , DN2YFE , DY2THP , 

#  DY2FEP,X12,X22,RK12,RK22) 

CALL  PC (DN 1XTH , DN 1XFE , DXl THF , DX 1 FEF , DN 1 YTH , DN 1 YFE , DY 1 THF , 

#  DY1 FEF, XI 1 , X21 , RKl 1 , RK21) 

CALL  PD (DX2THP , DX2FEP , DY2THP , DY2FEP , DZ2THP , DZ2FEP , XI 2 , 

//  PDX12,PDY12,PDZ12) 

CALL  PD (DX2THP , DX2FEP , DY2THP , DY2FEP , DZ2THP , DZ2FEP , X22 , 

#  PDX22 , PDY22 , PDZ22) 

CALL  PD(DX1THF,DX1FEF, DY1THF,DY1FEF,DZ1THF,DZ1FEF,X1 1 , 

#  PDX1 1 , PDY 1 1 , PDZ 11) 

CALL  PD(DX1THF,DX1FEF,DY1THF,DY1FEF,DZ1THF,DZ1FEF,X21 , 

#  PDX21 , PDY21 , PDZ21) 

CALL  UNIT (PDXl 1 , PDY1 1 , PDZ1 1 ,UNIT11,UX11,UY11,UZ11) 

CALL  UNIT (PDX2 1 , PDY2 1 , PDZ2 1 , UNIT21 , UX21 , UY21 , UZ21) 

CALL  UNIT (PDX 1 2 , PDY 1 2 , PDZ 12,UNIT12,UX12,UY12,UZ12) 

CALL  UNIT (PDX22 , PDY22 , PDZ22 , UNIT22 , UX22 , UY22 , UZ22) 

CALL  UTRAN 1 (UX 1 1 , U Y 1 1 , UZ 1 1 , UTX 1 1 , UTY 1 1 , UTZ 1 1 ) 

CALL  UTRAN 1 (UX21 ,UY21 , UZ21 , UTX21 , UTY21 , UTZ21) 

CALL  UTRAN2(UX12,UY12,UZ12,UTX12,UTY12,UTZ12) 

CALL  UTRAN2 (UX22 , UY22 , UZ22 , UTX22 , UTY22 , UTZ22) 

CALL  AXIS (SIGMDGjAXISA, AXISB , ALPH, ALPHDG) 

CALL  WRITE2 
FEP=FEP+DFEPDG:  CNST 
88  CONTINUE 
STOP 
END 

.  *****  SUBROUTINE  ANGLE  ***** 

SUBROUTINE  ANGLE (X) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/Al/CNST 

CNST2»2.DOO*CNST*180.DOO 

M-X/CNST2 

RM-M 

X=X-RM*CNST2 

RETURN 

END 

,  .  *****  SUBROUTINE  FCN1  ***** 

SUBROUTINE  FCN1 (X , F , N5 , PAR5) 


I 


I 
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IMPLICIT  REAL*8(A-H,0-Z) 

INTEGER  N5 

REAL*8  X(N5) ,F(N5) ,PAR5(N5) 

COMMON/Al/CNST 

COMMON/A2/DE 1 , DLI , RCP , RL , RCF, PHFDG 

COMMON/A3/SNPHF , CSPHF , SNR 1D1 , CSR 1D1 , SNBP , CSBP , SNR  1 D2 , CSR 1 D2 
COMMON /A5/SNGM1 , CSGMl , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 
COMMON /A6/ SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON /A7/SNDLTA.CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
C0MM0N/A12/UF,BF,RBP,UP,RMF1 ,RMP2 

COMMON /Al 8 /SNQP , CSQP , QPDG . SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON /A2 7 /XF 1 ,YF1 ,ZF1 ,XF2, YF2.ZF2 

COMMON/A34/FE1 , SNFE1 , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE22 1 , FE1 1 1 , 
ft  FEPDG,  FEFDG,  F111DG.F221DG 

COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
COMMON/ A37/RNXF2,RNYF2,RNZF2,RNXF1 ,RNYF1,RNZF1 
QF-QFO~X(3) 

SNQF-DSIN(QF) 

CSQF-DCOS(QF) 

QP-QPO-FEP 
SNQP-DS I N (QP) 

CSQP-DCOS(QP) 

FE1 1 1-X (4) 

FE221-X  (5) 

SN221-DSIN (FE221 ) 

CS221-DCOS (FE221) 

SNlll-DSIN(FElll) 

CSlll-DCOS(FElll) 

SNTF-DSIN(X(1) ) 

CSTF-DCOS(Xd)) 

SNTP-DS IN (X (2) ) 

CSTP-DCOS(X(2) ) 

THF-X(1)+QF 
THP-X (2) +QP 
THFDG- THF/CNST 
THPDG- THP/CNST 
SNTHF-DSIN (THF) 

CSTHF-DCOS (THF) 

SNTHP-DS I N (THP ) 

CSTHP-DCOS (THP) 

UFl-( (RCF*CSPHF/SNPHF-RL*SND1-DL1*SNR1D1/CSR1D1)*CSPHF*SNTF+ 
ft  BF*  (SNPHF*SNQF+CSPHF*SNTHF*  (RMF 1  -  SNR  1D1)/CSR1D1)-DE1*( 

ft  SNPHF-CSPHF*CSTF*SNR1D1/CSR1D1))/SNTF 

XH1- (CSGMl *CSR1D1 *CS1 1 1  "-SNGM1  *SNR  1D1 )  *  (RCF*CSPHF/SNPHF~UF1  * 
ft  CSPHF-RL*SND1)+BF* (CSGMl* (SNQF*SN1 1 1-CSQF*CS1 1 1* 

ft  SNR1D1 ) +SNGM1*CSQF*CSR 1D1 ) +DL1 * (SNGM1 *CSR 1D1 -CSGMl * 

ft  SNR 1D1 *CS 1 1 1 ) -SN 1 1 1 *DE1 *CSGM1 -UF 1 *SNPHF* (CSGMl * (SNTF* 

ft  SN111+CSTF*CS111*SNR1D1)-SNGM1*CSTF*CSR1D1) 

YH1-CSR1D1*SN111*(RCF*CSPHF/SNPHF-UF1*CSPHF-RL*SND1)-BF*( 
ft  SNQF*CS 1 1 1 +CSQF*SN 1 1 1 *SNR 1D1 ) -DL 1 *SNR 1D1 *SN 111 ♦DEI  * 
ft  CS1 1 1+UF1 *SNPHF* (SNTF*CS1 1 1~CSTF*SN 111 *SNR 1D1 ) 

ZH1A- (CSGMl *SNR1D1-SNGM1*CSR1D1*CS1  ID*  (RCF*CSPHF/SNPHF-UF1* 
ft  CSPHF) +RL*SND1*(SNGM1*CSR1D1*CS1 11 -CSGMl *SNR1D1)-BF*(SNGM1 

ft  * (SNQF*SN1 1 1~CSQF*CS 1 1 1 *SNR 1D1 ) -CSGMl *CSQF*CSR 1D1 ) ♦ 

ft  DL 1  * (SNGM 1 *SNR 1D1 *CS 1 1 1 +CSGM 1*CSR1D1)+DE1*SN111 *SNGM 1 
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ZH1»ZH1A+UF1*SNPHF*(SNGM1*(SNTF*SN1 1 1+CSTF*CS1 1 1 *SNR 1D1 ) +CSGM1 * 
it  CSTF’-CSRIDI) 

XF1C (CSDLTA*CSGM1**2+SNGM1**2) *XH1+SNDLTA*CSGM1,vYH1  + ( 1 . ODOO- 
//  CSDLTA)*CSGM1*SNGM1*ZH1 

YF1'-SNDLTA*CSGM1*XH1+CSDLTA*YH1+SNDLTA*SNGM1*ZH1 
ZF 1  -  ( 1 .  ODOO-CSDLTA)  *CSGM  1  *SNGM  1  *XH 1  -SNDLTA*SNGM1  * YH 1  + 
it  (CSDLTA*SNGMl**2-*-CSGMl**2)  *ZH1 

XF2a*RBP*  (CSPHP*SNTHP,'SND2/CSD2-SNPHP*SNQP)  *  (CSPHP*CSR2D2 
it  '"CSGM2,:CS22 1 /SNTP +CS PH P* SNR 2D2*SNGM2/ SNTP- SNR 2D2*SNPHP* 

it  CSGM2’:CS221,',CSTP/SNTP+SNPHP’',CSR2D2*SNGM2*CSTP/SNTP-SNPHP* 

it  CSGM2*SN22 1 )  +CSPHP*CSGM2*  (RCP’!CSPHP+RL*SNPHP*SND2)  *  (SNTP* 

it  SN22 1+CSTP*CS221*SNR2D2) 

XF2“XF2a+RBP:'CSGM2*  (SNQP*SN221+CSQP*CS221*SNR2D2) +RL*SND2* 
it  (CSCM2*CSR2D2*CS221+SNGM2*SNR2D2)-RBP*CSR2D2*CSQP*SNGM2 

it  +CSGM2*CSPHP*CSR2D2,tCS221,'ir  (RCP*SNPHP-RL*CSPHP*SND2)  -RCP* 

it  CSPHP*SNGM2* (CSPHP*CSR2D2*CSTP-SNPHP*SNR2D2) -RL*CSPHP*SND2 
it  *SNGM2*  (SNPHP*CSR2D2,tCSTP+CSPHP*SNR2D2) 

YF2=CSPHP::CSR2D2*SN22 1*  (RL*CSPHP*SND2-RCP*SNPHP)  +RBP*  (SNPHP* 
it  SNQP~CSPHP*SNTHP*SND2/CSD2)*(SNPHP*CS221-SNPHP*SNR2D2* 

it  SN22 1  *CSTP/SNTP+CSPHP*CSR2D2*SN221/SNTP)  -RBP*  (SNQP*CS221 

it  +CSOP*SN221  *SNR2D2)  -RL*SND2,rSN221*CSR2D2+CSPHP*  (RCP* 

it  CSPHP+RL*SNPHP*SND2)*(SNTP*CS221-CSTP*SNR2D2*SN221) 

ZF2A-RBP* (SNPHP*SNQP-CSPHP*SNTHP*SND2/CSD2) * (CSPHP*SNR2D2* 
it  CSGM2/SNTP-CSPHP*CSR2D2*SNGM2*CS221/SNTP+SNR2D2*SNPHP* 

it  SNGM2)'CS221*CSTP/SNTP+SNPHP*CSR2D2*CSGM2*CSTP/SNTP+SNPHP* 

it  SNGM2*SN221)-RBP*SNGM2*(SNQP*SN221-CSQP*CS221*SNR2D2) 

it  +RBP*C:R2D2*CSQP*CSGM2 

ZF2*ZF2a+SNGM2*CSPHP*  (SNTP*SN221+CSTP,rSNR2D2*CS221 )  *  (RCP*CSPHP 
it  +RL*SNPHP*SND2)+SNGM2*CSR2D2*CS221*(RCP*CSPHP*SNPHP-RL* 

it  SND2*CSPHP**2)  +RCP*CSGH2*CSPHP*  (CSPHP*CSR2D2*CSTP-SNPHP* 

it  SNR2D2)+RL*CSPHP*SND2*CSGM2*(CSPHP*SNR2D2+SNPHP*CSR2D2* 

it  CSTP)+RL*SND2*(SNGM2*CSR2D2*CS221-SNR2D2*CSGM2) 

RNXF2-CSPHP*CSGM2*SNTP*SN221+SNPHP*(SNGM2,tSNR2D2+CSGM2*CSR2D2* 
it  CS221)-CSPHP*CSTP*(SNGM2*CSR2D2-CSGM2*SNR2D2*CS221) 

RNYF2“CSPHP*CS221 *SNTP-SN221* (SNPHP *CSR2D2+CSPHP*SNR2D2*CSTP) 
RNZF2-CSPHP*CSTP*(CSGM2*CSR2D2+SNGK2*SNR2D2*CS221)+SNPHP*( 
it  SNGM2*CSR2D2*CS221-CSGM2*SNR2D2)+CSPHP*SNGM2*SN221*SNTP 

RNXH1  “-CSPHF*CSGM  1  *SNTF*SN1 1 1+SNPHF*  (SNGM1*SNR1D1+CSGM1* 
it  CSR 1D1 *CS 1 1 1 ) +CSPHF*CSTF*  (SNGM1  *CSR  1D1-CSGM1  *SNR  1D1  * 

it  CS 1 1 1 ) 

RNYH 1  “CSPHF*SNTF,!CS  1 1 1  +  SN1 1 1  * (SNPHF*CSR 1D1  ~CSPHF*SNR  1D1  *CSTF) 
RNZH1-CSPHF*CSTF*  (CSGM1  *CSR  1 D1 +  SNGM 1  *SNR  1D1  *CS  1 1 1 )  +SNPHF*  ( 
it  CSGM1 *SNR  1D1 -SNGM1 *CSR  1D1  *CS  1 1 1 )  +CSPHF*SNGMl*SNTF,:SNl  1 1 

RNXF1* (CSDLTA*CSGMl **2+SNGMl *<r2) *RNXH1 +SNDLTA*CSGM1 *RNYH 1 + 
it  ( 1 .  doo-csdlta)  *CSGM1  *SNGMl  *RNZH  1 

RNYF1«-SNDLTA*CSGM1*RNXH1+CSDLTA*RNYH1+SNDLTA*SNGM1*RNZH1 
RNZF1 - ( 1 . DOO-CSDLTA) *CSGN 1 *SNGM 1 *RNXH1 -SNDLTA*SNGM1 *RNYH 1 + 
it  (CSDLTA*SNGN1**2*CSGM1**2)  *RNZH1 

F(1)-XF1-XF2 
F(2)-YF1-YF2 
F  (3) “ZF1 -ZF2 
F(A)-RNXF1-RNXF2 
F(5)-RNYF1-RNYF2 
RETURN 
END 
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*****  SUBROUTINE  FCN2  ***** 

SUBROUTINE  FCN2 (X, F,N5, PAR5) 

IMPLICIT  REAL*8(A-H,0-Z) 

INTEGER  N5 

REAL*8  X(N5) ,F(N5) ,PAR5(N5) 

COMMON /Al/CNST 

COMMON /A2/ DEI , DL1 , RCP , RL, RCF, PHFDG 

C0MM0N/A3/SNPHF,CSPHF,SNR1D1,CSR1D1 ,SNBP,CSBP,SNR1D2,CSR1D2 
COMMON /A5/SNGM1 , CSGM 1 , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA , SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/Al 2/UF , BF , RBP , UP , RMF 1 , RMP2 

COMMON/ A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON /A27/XF1 , YF1 ,ZF1 , XF2, YF2 , ZF2 

COMMON/A3A/FE1 , SNFEI , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FE1 1 1 , 
it  FEPDG,  FEFDG,  F111DG,F221DG 

COMMON /A35/SNTHP , CSTHP , SNTHF , CSTHF, THPDG , THFDG 
QF-QFO+XC3) 

SNQF*DSIN (QF) 

CSQF-DCOS(QF) 

QP=QPO+FEP 
SNQP=DSIN (QP) 

CSQP=DCOS (QP) 

FE1 1 1*X (4) 

FE221“X  (5) 

SN221“DSIN (FE221) 

CS221*DCOS (FE221) 

SNlll-DSIN(FEIII) 

CSlll-DCOS(FElll) 

SNTF-DSIN(X(1)) 

CSTF-DCOS (X(l)) 

SNTP”DSIN (X (2) ) 

CSTP-DCOS (X (2) ) 

THF“X(1)-QF 
THP'X (2) -QP 
THFDG-THF/CNST 
THPDG-THP/CNST 
SNTHF-DSIN(THF) 

CSTHF-DCOS(THF) 

SNTHP-DSIN(THP) 

CSTHP-DCOS(THP) 

UF1-((RCF*CSPHF/SNPHF~RL*SND1-DL1*SNR1D1/CSR1D1)*CSPHF*SNTF+ 
H  BF*(-SNPHF*SNQF+CSPHF*SNTHF*(RMF1-SNR1D1)/CSR1D1)-DE1*( 

it  SNPHF-CSPHF*CSTF*SNR  1D1  /CSR 1D1 ) )  /SNTF 
XH1«(CSGM1*CSR1D1*CS111+SNGM1*SNR1D1)*(RCF*CSPHF/SNPHF-UF1* 
it  CSPHF-RL*SND1)  +BF*  (CSGM1  *  (~SNQF*SNl  1 1~CSQF*CS1 1 1  * 

it  SNR1D1)  +SNGMl*CSQF*CSR  1D1 )  +DL1  *  (SNGM1  *CSR  1DI-CSGM1  * 

it  SNR1D1*CS1 1 1)  -SN11 1*DE1  *CSGM1-UF1*SNPHF*  (CSGMl*  (SNTF* 
it  SN1 1 1+CSTF*CS1 1 1*SNR  1D1)  -SNGM1  *CSTF*CSR  1D1) 

YH1-CSR1D1*SN1 1 1* (RCF*CSPHF/SNPHF-UF1*CSPHF-RL*SND1) -BF* ( 
if  _SNQF*CS1 1  l+CSQF*SNl  1 1  *SNR  1D1)  ~DL1*SNR  1D1*SN1 1 1+DE1* 
it  CS1 1 1+UF1*SNPHF*  (SNTF*CS1 1 1-CSTF*SN1 1 1*SNR1D1) 

ZHlA- (CSGMI*SNR1D1-SNGM1*CSR1D1*CS111)*(RCF*CSPHF/SNPHF-UF1* 
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//  CSPHF)  +RL*SND1 *  (SNGM1 *CSR 1D1  *CS 1 1 1 -CSGM1 *SNR 1D1 ) -BF* (SNGM1 

it  *(-SNQF*SNlll-CSQF*CSlll*SNRlDl)-CSGMl*CSQF*CSRlDl)+ 

if  DL1*  (SNGM1*SNR1D1*CS1 1 1+CSGM1*CSR1D1)  +DE1*SN1  1  l*SNGMl 

ZH 1=ZH1 A+UF1 *SNPHF*  (SNGM1*  (SNTF*SN1 1 1+CSTF*CS 1 1 1  *SNR 1D1 )  +CSGM1  * 
if  CSTF*CSR1D1) 

XF1=  (CSDLTA*CSGM1 **2+SNGMl **2)  *XH1+SNDLTA*CSGM1 *YH 1+  ( 1 . 0D00- 
//  CSDLTA)  *CSGM1  *SNGM1*ZH1 

YF1=-SNDLTA*CSGM1*XH1+CSDLTA*YH1+SNDLTA*SNGM1*ZH1 
2  F 1 = ( 1 . ODOO-CSDLTA) *CSGM 1 *SNGM1 *XH 1 -SNDLTA" SNGM 1  * YH 1  + 
if  (CSDLTA*SNGM1**2+CSGM1**2)*2H1 

XF2A-RBP*  (CSPHP*SNTHP*SND2/CSD2+SNPHP*SNQP)  *  (CSPHP*CSR2D2 
if  *CSGM2*CS221/SNTP+CSPHP*SNR2D2*SNGM2/  SNTP- SNR 2D2*SNPHP* 

it  CSGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*SNGM2*CSTP/SNTP-SNPHP* 

if  CSGM2*SN221)+CSPHP*CSGM2*(RCP*CSPHP+RL*SNPHP*SND2)  *  (SNTP* 

if  SN22 1+CSTP*CS22 1  *SNR2D2) 

XF2=XF2A+RBP*CSGM2*(-SNQP*SN221+CSQP*CS221*SNR2D2)+RL*SND2* 
it  (CSGM2*CSR2D2*CS221+SNGM2*SNR2D2)-RBP*CSR2D2*CSQP*SNGM2 

it  +CSGM2*CSPHP*CSR2D2*CS221*  (RCP*SNPHP-RL*CSPHP*SND2) -RCP* 

it  CSPHP*SNGM2* (CSPHP*CSR2D2*CSTP-SNPHP*SNR2D2) -RL*CSPHP*SND2 

it  *SNGM2*  (SNPHP*CSR2D2*CSTP+CSPHP*SNR2D2) 

YF2=CSPHP*CSR2D2*SN221*(RL*CSPHP*SND2-RCP*SNPHP)+RBP*(-SNPHP* 
it  SNQP-CSPHP*SNTHP*SND2/CSD2)  *  (SNPHP*CS221-SNPHP*SNR2D2* 

it  SN221*CSTP/SNTP+CSPHP*CSR2D2*SN221/SNTP)-RBP*(-SNQP*CS221 

it  +CSQP*SN221*SNR2D2)-RL*SND2*SN221*CSR2D2+CSPHP*(RCP* 

it  CSPHP+RL*SNPHP*SND2)*(SNTP*CS221-CSTP*SNR2D2*SN221) 

2F2A=RBP* (-SNPHP*SNQP-CSPHP*SNTHP*SND2/CSD2) * (CSPHP*SNR2D2* 
it  CSGM2/SNTP-CSPHP*CSR2D2*SNGM2*CS221/SNTP+SNR2D2*SNPHP* 

it  SNGM2*CS221*CSTP/SNTP+SNPHP*CSR2D2*CSGM2*CSTP/SNTP+SNPHP* 

it  SNGM2*SN221)-RBP*SNGM2*(-SNQP*SN221-CSQP*CS221*SNR2D2) 

it  +RBP*CSR2D2*CSQP*CSGM2 

ZF2»ZF2A+SNGM2*CSPHP* (SNTP*SN221+CSTP*SNR2D2*CS221) * (RCP*CSPHP 
it  +RL*SNPHP*SND2)  +SNGM2*CSR2D2*CS221 *  (RCP*CSPHP*SNPHP-RL* 

it  SND2*CSPHP**2) +RCP*CSGM2*CSPHP* (CSPHP*CSR2D2*CSTP~SNPHP* 

it  SNR2D2)  +RL*CSPHP*SND2*CSGM2*  (CSPHP*SNR2D2+SNPHP*CSR2D2* 

it  CSTP)+RL*SND2*(SNGM2*CSR2D2*CS221-SNR2D2*CSGM2) 

RNXF2-CSPHP*CSGM2*SNTP*SN221+SNPHP*(SNGM2*SNR2D2+CSGM2*CSR2D2* 
if  CS221)  ~CSPHP*CSTP*  (SNGM2*CSR2D2-CSGM2*SNR2D2*CS221) 

RNYF2=CSPHP*CS221*SNTP-SN221*(SNPHP*CSR2D2+CSPHP*SNR2D2*CSTP) 
RNZF2-CSPHP*CSTP* (CSGM2*CSR2D2+SNGM2*SNR2D2*CS221) +SNPHP* ( 
it  SNGM2*CSR2D2*CS221-CSGM2*SNR2D2)+CSPHP*SNGM2*SN221*SNTP 

RNXH1— CSPHF*CSGM1*SNTF*SN111+SNPHF*(SNGM1*SNR1D1+CSGM1* 
it  CSR1D1  *CS  1 1 1 )  +CSPHF*CSTF*  (SNGM1 *CSR 1D1-CSGM1 *SNR 1D1  * 

it  CS111) 

RNYH1-CSPHF*SNTF*CS1 1 1+SN1 1 1* (SNPHF*CSR1D1-CSPHF*SNR1D1*CSTF) 
RNZH1*CSPHF*CSTF* (CSGM1 *CSR 1D1+SNGM1 *SNR1D1 *CS1 1 1) +SNPHF* ( 
it  CSGM1*SNR1D1-SNGM1*CSR1D1*CS111)+CSPHF*SNGM1*SNTFaSN111 

RNXF1-(CSDLTA*CSGM1**2+SNGM1**2)*RNXH1+SNDLTA*CSGM1*RNYH1+ 
it  ( 1 .  DOO-CSDLTA)  *CSGMl  *SNGM1  *RNZHl 

RNYF1-~SNDLTA*CSGM1*RNXH1+CSDLTA*RNYH1+SNDLTA*SNGM1*RNZH1 
RNZF1- ( 1 . DOO-CSDLTA) *CSGMl*SNGMl*RNXHl-SNDLTA*SNGMl*RNYH 1 + 
it  (CSDLTA*SNGM1**2+CSGM1**2)*RNZH1 

F (1) “XF1-XF2 
F  (2) “YF1-YF2 
F  (3) “ZF1-ZF2 
F  (A) “RNXF1-RNXF2 
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non  non  non 


F(5)'RNYF1-RNYF2 

RETURN 

END 

*****  SUBROUTINE  FNORM  ***** 

SUBROUTINE  FNORM (DN 1XTH , DN 1XFE , DN 1 YTH , DN 1 YFE , DN 1 ZTH , DN 1ZFE) 
IMPLICIT  REAL*8(A-H,0-Z) 

C0MM0N/A3/SNPHF,CSPHF,SNR1D1,CSR1D1,SNBP,CSBP,SNR1D2,CSR1D2 
COMMON /a6/ SNTP , CSTP , TP , TPDG , SNTF, CSTF , TF , TFDG 
C0MM0N/A12/UF,BF,RBP,UP,RMF1,RMP2 

COMMON /a3 A /FE1 , SNFE1 , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FE1 1 1 , 
it  FEPDG,  FEFDG,  F111DG.F221DG 

DN1XTH«-CSPHF*SNFE1*CSTF+CSPHF*CSFE1*SNR1D1*SNTF 
DN1XFE— SNPHF*SNFE1*CSR1D1/RMF1-CSPHF*(CSFE1*SNTF/RMF1+SNFE1* 
it  CSTF)-CSPHF*SNR1D1*(-SNFE1*CSTF/RMF1-CSFE1*SNTF) 

DN 1 YTH=CSPHF*CSFE1 *CSTF+CSPHF*SNFE1*SNR1D1 *SNTF 
DN  1  YFE“SNPHF*CSFE1 *CSR 1D1/RMF1+CSPHF*  (_SNFE1*SNTF/RMF1+CSFE1  * 
it  CSTF)-CSPHF*SNR1D1*(CSFE1*CSTF/RMF1-SNFE1*SNTF) 

DN1ZTH--CSPHF*CSR1D1*SNTF 
DN1ZFE“-CSPHF*CSR1D1*SNTF 
RETURN 
END 

*****  SUBROUTINE  PNORM  ***** 

SUBROUT I NE  PNORM (DN2XTH , DN2XFE , DN2 YTH , DN2 YFE , DN2ZTH , DN2ZFE) 
IMPLICIT  REAL*8(A-H,0-Z) 

C0MM0N/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA, CSDLTA, SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/ A12/UF , BF , RBP , UP , RMF1 , RMP2 

COMMON/ A3A/ FE 1 , SNFE 1 , CSFEl , FE2 , SNFE2 , CSFE2 , FEP , FE22 1 , FE 1 1 1 , 
it  FEPDG,  FEFDG,  FI  1  IDG,  F221DG 

DN2XTH«CSPHP*(SNFE2*CSTP-CSFE2*SNR2D2*SNTP) 

DN2XFE— SNFE2*CSR2D2*SNPHP/RMP2+CSPHP*(CSFE2*SNTP/RMP2+SNFE2* 
it  CSTP)  +CSPHP*SNR2D2*  (-SNFE2*CSTP/RMP2~CSFE2*SNTP) 

DN2YTH-CSFE2*CSPHP*CSTP+SNFE2*SNR2D2*CSPHP*SNTP 
DN2YFE— CSFE2*CSR2D2*SNPHP/RMP2+CSPHP*(-SNFE2*SNTP/RMP2+ 
it  CSFE2*CSTP)  -CSPHP*SNR2D2*  (CSFE2*CSTP/RMP2-SNFE2*SNTP) 

DN2ZTH— CSPHP*CSR2D2*SNTP 
DN2ZFE— CSPHP*CSR2D2*SNTP 
RETURN 
END 

*****  SUBROUTINE  DUP  ***** 

SUBROUTINE  DUP (UPP , DUPTHP , DUPFEP) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/ A2/DE1,DL1,RCP,RL,RCF,PHFDG 

COMMON/ A5 / SNGM 1 , CSGM 1 , GMA 1 , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 
COMMON/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA , CSDLTA , SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/Al 2/UF , BF , RBP , UP , RMF 1 , RMP2 

COMMON/ A18 /SNQP , CSQP , QPDG, SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/ A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
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•f1  >■'  >i 


non  non 


UPP=RBP* (SNPHP*SNQP~CSPHP*SNTHP*SND2/CSD2) /SNTP+CSPHP* (RCP* 
#  CSPHP/SNPHP+RL*SND2) 

DUPTHP-RBP*SNQP*(~SNPHP*CSTP+CSPHP*SND2/CSD2)/SNTP**2 
DUPFEP*-RBP*SNTHP*  (SNPHP~CSPHP*CSTP*SND2/CSD2)  /SNTP**2 
RETURN 
END 


*****  SUBROUTINE  DUF  ***** 

SUBROUTINE  DUF (UFF, DUFTHF , DUFFEF) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON / A2 /DE 1 , DL1 , RCP , RL , RCF , PHFDG 

C0MM0N/A3/SNPHF,CSPHF,SNR1D1,CSR1D1,SNBP,CSBP,SNR1D2,CSR1D2 
COMMON /A5/SNGM1 , CSGM1 ,GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSDl 
C0MM0N/a6/SNTP,CSTP,TP,TPDG,SNTF,CSTF,TF,TFDG 
COMMON/A12/UF , BF , RBP , UP , RMF1 , RMP2 

COMMON/Al 8/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG , THFDG 
UFF=  (BF* (SNPHF*SNQF+CSPHF*SNTHF* (RMF1 /CSR 1D1 -SNR 1D1 /CSR 1D1 ) ) - 
//  DEI  *  (SNPHF~CSPHF*CSTF*SNR  ID1  /CSR  1D1 )  +CSPHF*SNTF*  (RCF* 

if  CSPHF/SNPHF-RL*SND1-DL1*SNR1D1/CSR1D1))/SNTF 
DUFTHF=(-BF*SNQF*(SNPHF*CSTF+CSPHF*(RMF1/CSR1D1-SNR1D1/CSR1D1 
if  ))-DEl*(CSPHF*SNRlDl/CSRlDl-SNPHF*CSTF))  /SNTF**2 

DUFFEF= (SNPHF+CSPHF*CSTF* (RMF1 /CSR 1D1-SNR 1D1 /CSR1D1))*(_BF* 
if  SNTHF/SNTF**2)  +DE1  *  (SNPHF*CSTF~CSPHF*SNR  1D1  /CSR1D1)  / 

if  SNTF**2 

RETURN 
END 

*****  SUBROUTINE  DR2  ***** 

SUBROUTINE  DR2 (UPP , DUPTHP , DUPFEP , DX2THP , DX2FEP . DY2THP , DY2FEP , 
if  DZ2THP ,  DZ2FEP) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/A2/DE1 ,DL1 , RCP ,RL, RCF, PHFDG 

COMMON/A5/SNGM1 , CSGMI , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSDl 
COMMON/ A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A7/SNDLTA , CSDLTA , SNPHP , CSPHP , SNR2D2 , CSR2D2 , DLTADG 
COMMON/Al 2/UF , BF , RBP , UP , RMF1 , RMP2 

COMMON/A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/ A34/FE1 , SNFE1 , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FE1 1 1 , 
if  FEPDG,  FEFDG,  FI  1  IDG,  F221DG 

DX2THP— CSPHP*CSFE2*CSR2D2*DUPTHP+SNPHP*SNFE2*(UPP*CSTP+SNTP 
if  *DUPTHP)  +CSFE2*SNR2D2*SNPHP*  (~UPP*SNTP+CSTP*DUPTHP) 

DX2FEP— RCP*CSR2D2*SNFE2* (CSPHP/SNPHP) /RMP2~CSPHP*CSR2D2* ( 
if  -UPP*SNFE2/RMP2+CSFE2*DUPFEP)  -RL*SND2*CSR2D2*SNFE2/RMP2 

if  -RBP*  (-SNFE2*CSQP  +  SNQP*CSFE2/RMP2)  +RBP*SNR2D2*  (CSFE2* 

if  SNQP-CSQP*SNFE2/RMP2) *SNPHP* (UPP* (SNFE2*CSTP+SNTP*CSFE2 

if  /RMP2)  +SNTP*SNFE2*DUPFEP) 

DX2FEP-DX2FEP+SNPHP*SNR2D2* (UPP* (-CSFE2*SNTP-CSTP*SNFE2/RMP2) + 
if  CSFE2*CSTP*DUPFEP) 

DY2THP«SNFE2*CSR2D2*CSPHP*DUPTHP+SNPHP*CSFE2*(SNTP*DUPTHP+UPP 
if  *CSTP)  -SNPHP*SNFE2*SNR2D2*  (-UPP*SNTP+CSTP*DUPTHP) 

DY2FEP— RCP*CSR2D2*CSFE2* (CSPHP/SNPHP) /RMP2+CSPHP*CSR2D2* (UPP* 
if  CSFE2/RMP2+SNFE2*DUPFEP)  -RL*SND2*CSR2D2*CSFE2/RMP2-RBP* 
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//  (-CSFE2*CSQP-SNFE2*SNQP/RMP2)  -RBP*SNR2D2*  (SNFE2*SNQP+ 

//  CSQP*CSFE2/RMP2)  +SNPHP*  (UPP*  (CSFE2*CSTP-SNFE2*SNTP/RMP2 

#  )+CSFE2*SNTP*DUPFEP) 

DY2FEP=DY2FEP-SNPHP*SNR2D2*(UPP*(-SNFE2*SNTP+CSTP*CSFE2/RMP2) 

//  +SNFE2*CSTP*DUPFEP) 

DZ2THP=SNR2D2*CSPHP*DUPTHP+SNPHP*CSR2D2*(-UPP*SNTP+CSTP* 

#  DUPTHP) 

DZ2FEP=CSPHP*SNR2D2*DUPFEP+RBP*CSR2D2*SNQP+SNPHP*CSR2D2*(-UPP 
//  *SNTP+CSTP*DUPFEP) 

RETURN 

END 

C.  .  . 

C...  *****  SUBROUTINE  DR  1  ***** 

C.  .  . 

SUBROUTINE  DRl (UFF,DUFTHF,DUFFEF,DX1THF,DX1FEF,DY1THF,DY1FEF, 

//  DZ1THF.DZ1FEF) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON / A2 / DE 1 , DLl , RCP , RL , RCF , PHFDG 

COMMON/ A3/ SNPHF , CSPHF , SNR1D1 , CSR 1D1 , SNBP , CSBP , SNR1D2 , CSR1D2 
COMMON/A5/SNGM1 , CSGM 1 , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 
COMMON/A6/ SNTP , CSTP , TP , TPDG , SNTF, CSTF , TF , TFDG 
COMMON/A12/UF, BF , RBP , UP , RMF1 , RMP2 

COMMON/A18/SNQP , CSQP , QPDG , SNQF , CSQF , QFDG , QP , QF , QPO , QFO 
COMMON/A34/FE1 , SNFE1 , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE22 1 , FE1 1 1 , 

//  FEPDG,  FEFDG,  FI  1  IDG,  F221DG 

DX1THF=-CSPHF*CSFE1*CSR1D1*DUFTHF-SNPHF*SNFE1*(UFF*CSTF+SNTF* 

//  DUFTHF)  -SNPHF*CSFE1*SNR1D1*  (-UFF*SNTF+CSTF*DUFTHF) 

DX1FEF=-RCF*CSR1D1*SNFE1* (CSPHF /SNPHF) /RMF1~CSPHF*CSR1D1* ( 

//  -UFF*SNFE1/RMF1+CSFE1*DUFFEF)+RL*SND1*CSR1D1*SNFE1/RMF1 

#  +BF* (-SNFE 1 *CSQF+SNQF*CSFE1 /RMF1 ) -BF*SNR 1 Dl * (CSFE1 *SNQF 

#  -CSQF*SNFE1/RMF1)+DL1*SNFE1*SNR1D1/RMF1-DE1*CSFE1/RMF1 
DX1 FEF=DXlFEF~SNPHF* (UFF* (SNFE1*CSTF+SNTF*CSFE1/RMF1) + 

#  SNFE 1 *SNTF*DUFFEF) -SNPHF*SNR 1D1 * (UFF* (-CSFE 1 *SNTF~CSTF* 

#  SNFEl/RMFl) +CSFEl*CSTF*DUFFEF) 
DY1THF=-CSPHF*SNFE1*CSR1D1*DUFTHF+SNPHF*CSFE1*(UFF*CSTF+SNTF* 

#  DUFTHF) -SNPHF*SNFE1*SNR1D1*(-UFF*SNTF+CSTF*DUFTHF) 
DY1FEF=RCF*CSR1D1*CSFE1*(CSPHF/SNPHF)/RMF1-CSPHF*CSR1D1*( 

//  UFF*CSFE1/RMF1+SNFE1*DUFFEF)-RL*SND1*CSR1D1*CSFE1/RMF1- 

#  BF*(-CSFE1*CSQF-SNFE1*SNQF/RMF1)-BF*SNR1D1*(SNFE1*SNQF 

#  +CSFE1*CSQF/RMF1)  -DL1*CSFE1*SNR1D1/RMF''  DF 1  *SNFE1  /RMF 1 
DY1FEF”DY1FEF+SNPHF* (UFF* (CSFE1*CSTF-SNFE1 *SNTF/RMF1) +CSFE1* 

#  SNTF*DUFFEF)-SNPHF*SNR1D1*(UFF*(-SNFE1*SNTF+CSFE1*CSTF/ 

#  RMF1) +SNFE1*CSTF*DUFFEF) 

DZ1THF— CSPHF*SNR1D1*DUFTHF+SNPHF*CSR1D1*(-UFF*SNTF+CSTF* 

#  DUFTHF) 

DZ 1 FEF— CSPHF*SNR 1D1 *DUFFEF+BF*CSR 1D1 *SNQF+SNPHF*CSR 1D1 * ( 

#  -UFF*SNTF+CSTF*DUFFEF) 

RETURN 

END 

C.  . . 

*****  SUBROUTINE  PRINCIPAL  CURVATURE  ***** 

C. .  • 

SUBROUTINE  PC (Dl , D2 , D3 , DA , D5 , D6 , D7 , D8 , XI , X2 , RK1 ,RK2) 

IMPLICIT  REAL*8(A-H,0-Z) 

A-D1*D7-D5*D3 
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non  non  non  non  non 


B»D1*D8+D2*D7-D5*D4-D6*D3 

C=D2*D8-D6*D4 

Xl= (-B-DSQRT (B**2~4 . 0D00*A*C) ) / (2 . 0D00*A) 
X2=  (-B+DSQRT  (B**2-4 . 0D00*A*C) )  /  (2 . 0D00*A) 
RK1=- (Dl*Xl+D2) / (D3*X1+D4) 

RK2=- (D1*X2+D2) / (D3*X2+D4) 

RETURN 

END 


*****  SUBROUTINE  PRINCIPAL  DIRECTION  ***** 

SUBROUTINE  PD (Al , A2 , A3 , A4 , A5 , A6 , PCO , PD1 , PD2 , PD3) 
IMPLICIT  REAL*8(A-H,0-Z) 

PDl=Al*PC0+A2 

PD2=A3*PC0+A4 

PD3*A5*PC0+A6 

RETURN 

END 


*****  SUBROUTINE  UNIT  VECTOR  OF  PRINCIPAL  DIRECTION  ***** 

SUBROUTINE  UNIT (El , E2 , E3 , UNITO , UNIT1 ,UNIT2,UNIT3) 

IMPLICIT  REAL*8(A-H,0-Z) 

UNITO=DSQRT (El**2+E2**2+E3**2) 

UNI T1=E1 /UNITO 

UNIT2=E2/UNIT0 

UNIT3=E3/UNIT0 

RETURN 

END 


*****  SUBROUTINE  UTRAN1  ***** 

SUBROUTINE  UTRAN1 (B1 , B2 , B3 , UTl 1 , UT21 , UT31) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/A5/SNGM1 , CSGM1 , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 

UTl 1=B1*CSGM1+B3*SNGM1 

UT21=B2 

UT31=B1* (-SNGM1) +B3*CSGMl 

RETURN 

END 

*****  SUBROUTINE  UTRAN2  ***** 

SUBROUTINE  UTRAN2 (C 1 , C2 , C3 , UT 1 2 , UT22 , UT32) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/A5/SNGM1 , CSGM1 , GMAl , SNGM2 , CSGM2 , SND2 , CSD2 , SND1 , CSD1 

UT12-C1*CSGM2-C3*SNGM2 

UT22-C2 

UT32-C 1 *SNGM2+C3*CSGM2 

RETURN 

END 

*****  SUBROUTINE  ELLIPSE  AXIS  ***** 

SUBROUTINE  AXIS(SIGMDG, AXISA, AXISB,ALPH,ALPHDG) 
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non 


IMPLICIT  REAL*8(A-H,0-Z) 

REAL *8  SIGM(4) 

COMMON /Al /CNST 

COMMON/ A36/RK12 , RK22 ,RK11,RK21,UTX11,UTX12,UTY11,UTY12, UTZl 1 , 

#  UTZ 12 , UTX21 , UTX22 , UTY21 , UTY22 , UTZ21 , UTZ22 
COMMON / A3 7 /RNXF2 ,RNYF2,RNZF2 ,RNXFl ,RNYFl ,RNZFl 
COMMON/A38/DEF , S IGMGG, ALPHGG , AXIAA , AXIBB , ERROR 
AK2=RK12+RK22 

AK1=RK1 1+RK21 
GI=RKl 1-RK21 
G2=RK12-RK22 

C...  FOR  LEFT-HAND  GEAR 

CS21=(UTX11*UTX12+UTY11*UTY12+UTZ11*UTZ12) 

IF(lI.EQ.2)GO  TO  1208 

SN21=UTXll*RNYF2*UTZ12+UTX12*RNZF2*UTYll+RNXF2*UTYl2*UTZll~ 

//  (UTX 1 2*RN Y  F2  *UTZ  1 1  +RNXF2*UT  Y 1 1  *UTZ  1 2+UTX 1 1  *RNZF2*UTY  1 2) 

GO  TO  1209 

1208  SN21=-(UTX11*RNYF2*UTZ12+UTX12*RNZF2*UTY11+RNXF2*UTY12*UTZ11) 

//  +(UTX12*RNYF2*UTZ11+RNXF2*UTY11*UTZ12+UTX11*RNZF2*UTY12) 

1209  SIGM21=DATAN2 (SN21 ,CS21) 

5  CALL  ANGLE (SIGM21) 

CSSIGM=DCOS (SIGM21) 

SNSIGM=DSIN (SIGM21) 

SIGMDG=SIGM21/CNST 

AA=(AKl-AK2~DSQRT(Gl**2-2.0D00*Gl*G2*DCOS(2.D00*SIGM21)+G2**2) 
//  )/4.0D00 

BB= (AK1-AK2+DSQRT (Gl**2-2 . D00*G1 *G2*DCOS (2.D00*SIGM21)+G2**2) ) 

#  /4.0D00 

AXISA=DSQRT (DABS (DEF/AA) ) 

AXISB=DSQRT (DABS (DEF/BB)  ) 

RATIO=AXISA/AXISB 
FF“G1*DSIN (2.0D00*SIGM21) 

HH=G2-G1*DC0S(2.D00*SIGM21) 

ALPH2=DATAN2 (FF, HH) 

ALPH=ALPH2/2 . 00D00 
CALL  ANGLE (ALPH) 

ALPHDG=ALPH/CNST 
RETURN 
END 

*****  SUBROUTINE  WRITE2  ***** 

SUBROUTINE  WRITE2 
IMPLICIT  REAL*8(A-H,0-Z) 

COMMON/A6/SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON /A2 7 /XFl , YF1 , ZFl , XF2 , YF2 , ZF2 

COMMON/ A34/FE 1 , SNFE 1 , CSFE1 , FE2 , SNFE2 , CSFE2 , FEP , FE221 , FEl 1 1 , 

#  FEPDG, FEFDG, FI  1  IDG , F221DG 
COMMON/A35/SNTHP , CSTHP , SNTHF , CSTHF , THPDG, THFDG 
COMMON/A36/RK 1 2 , RK22 , RK1 1 , RK2 1 , UTX 1 1 , UTX1 2 , UTY 1 1 , UTY 1 2 , UTZ 1 1 , 

#  UTZ 1 2 , UTX2 1 , UTX22 , UTY2 1 , UTY22 , UTZ2 1 , UTZ22 
COMMON/A38/DEF , SIGMDG , ALPHDG , AXISA , AXI SB , ERROR 
WRITE (6, 70002) 

70002  FORMAT (1H  , ' **RESULT  OF  KINEMATIC  ERROR  &  BEARING  CONTACT**’/) 
WRITE (6, 70003) FEPDG, FEFDG, THPDG, THFDG, TPDG, TFDG, F221DG, FI 1  IDG 


70003  FORMAT (1H  FEPDG= ' ,G20. 12, 10X, ' FEFDG= ' ,G20. 12, ' 

#  1H  , '  THPDG= ' , G20 . 1 2 , 1 OX , ' THFDG= ' , G20 .12,/ 

//  1H  TPDG  =' ,G20. 12, 10X, 'TFDG  =’,G20.12,/ 

#  1H  F221DG= ' ,G20. 12 ,9X, 1  FI 1 1DG= ' , G20. 12) 

WRITE (6 , 70004) XF1 , YF1 , ZF1 , XF2 , YF2 , 2F2 

70004  FORMAT (1H  , 'XF1=' ,G20. 12, 3X, ' YFl= ' ,G20. 12 , 3X, ’ZF1=’ ,G20. 12,/ 

#  1H  , ’XF2=' .G20.12.3X, *YF2=' ,G20.12,3X, 'ZF2=' .G20.12/) 
WRITE (6, 70005) ERROR 

70005  FORMAT (1H  , ' KINEMATIC  ERROR= ' ,G20. 12/) 

WRITE(6, 70006) 

70006  FORMAT (1H  ,'RESULT  OF  PRINCIPAL  CURVATURE') 

WRITE (6 , 70007) RK12 , RK22 , RK1 1 , RK2 1 

70007  FORMAT (1H  , 'GEAR  :  RK12= ' , G20. 12, 3X, ' RK22= ' , G20. 12, / 

//  1H  , ’PINION:  RK1 1=  '  , G20 . 12, 3X,  ' RK21=  '  , G20. 12/) 

WRITE (6 , 70008) UTX1 1 , UT Y 1 1 , UTZ 1 1 , UTX2 1 , UTY2 1 , UTZ2 1 

70008  FORMAT (1H  , 'PINION  UNIT  VECTOR  OF  PRINCIPAL  DIRECTION  ',/ 

#  1H  , 'UTX11=' ,G20.12, 'UTY11-' ,G20. 12, 'UTZ11=' ,G20. 12,/ 

#  1H  , 'UTX21=',G20.12,,UTY21=',G20.12, ’UTZ21=’,G20.12/) 

WR I TE (6 , 70009) UTX 1 2 , UTY 1 2 , UTZ 1 2 , UTX22 , UTY22 , UTZ22 

70009  FORMAT ( 1H  , ' GEAR  UNIT  VECTOR  OF  PRINCIPAL  DIRECTION’,/ 

#  1H  , ' UTX 12= ' , G20. 12 , ’ UTY12= ' ,G20. 12, ' UTZ 12= ' , G20. 12 , / 

//  1H  , 'UTX22=' ,G20.12, 'UTY22=' ,G20. 12, ’UTZ22=’ ,G20. 12/) 

WRITE (6, 70010) SIGMDG,ALPHDG,AXISA,AXISB 

70010  FORMAT (1H  .'DIRECTION  L  DIMENSION  OF  ELLIPSE’,/ 

#  1H  , ' SIGMDG= ' ,G20. 12, 10X, ' ALPHDG= ' , G20. 12 , / 

#  1H  , 'AXISA=' ,G20.12,11X, *AXISB=' ,G20.12,/ 

#  1H  ,  72 ( ' . ' ) /) 

RETURN 

END 

C.  .  . 

C...  *****  SUBROUTINE  HEADl  ***** 

C.  .  . 

SUBROUTINE  HEADl(JJ.Il) 

IMPLICIT  REAL*8(A-H,0-Z) 

WRITE (6, 60000) 

WRITE(6, 60001) 

WRITE (6,60002) 

IF ( J J  .EQ.  1)THEN 
WRITE(6, 60005) 

ELSE 

WRITE(6, 60006) 

END  IF 

WRITE(6, 60002) 

IF(II.EQ.2)  GO  TO  1 
WRITE (6, 60003) 

GO  TO  2 

1  WRITE (6, 60004) 

2  WRITE (6, 60002) 

WRITE(6, 60001) 

WRITE (6, 60000) 


60000  FORMAT (1H1 , '  ') 

60001  FORMAT (1H  '  ************************************************ ' ) 

60002  FORMAT (1H  , 

*’) 

60003  FORMAT (1H 

CONVEX 

PART  ANALYSIS 

*') 

60004  FORMAT (1H 

CONCAVE 

PART  ANALYSIS 

*') 

119 


u  u  o  u 


60005  FORMAT (1H  LEFT-HAND  GEAR  *') 

60006  FORMAT (1H  RIGHT-HAND  GEAR  *') 

RETURN 

END 


. ..  *****  WRITEO(II)  ***** 

SUBROUTINE  WRITEO(II) 

IMPLICIT  REAL*8(A-H,0-2) 

C0MM0N/A2/DE1 , DL1 , RCP , RL , RCF , PHFDG 
COMMON /A6/ SNTP , CSTP , TP , TPDG , SNTF , CSTF , TF , TFDG 
COMMON/A33/PHP1 , PHF1 , PHP2 , PHF2 , PHP1DG , PHF1DG, PHP2DG , PHF2DG 
IF(II.EQ.2)G0  TO  1 
PP-PHP1DG 
PF-PHF1DG 
GO  TO  90000 
1  PP“PHP2DG 
PF=PHF2DG 

90000  WRITE (6,90001) 

90001  FORMAT (1H  , 'GEAR  PARAMETER 27X P INION  PARAMETER’) 

WRITE (6 , 90002) PP , PF , RCP , RCF , TPDG , TFDG 

90002  FORMAT (1H  , 'PHPDG= ' ,G20. 12, 15X, 'PHFDG-' ,G20. 12,/ 

//  1H  ,'RCP  “  ’  ,G20. 12, 15X,  '  RCF  -'.G20.12,/ 

#  1H  , 'TPDG  -' ,G20.12,15X, 'TFDG  .G20.12) 

RETURN 

END 
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A  new  method  for  generation  of  spiral  bevel  gears  Is  proposed.  The  main  features 
of  this  method  are  as  follows:  (1)  the  gear  tooth  surfaces  are  conjugated  and  can 
transform  rotation  with  zero  transmission  errors;  (11)  the  tooth  bearing  contact 
Is  localized;  (111)  the  center  of  the  Instantaneous  contact  ellipse  moves  In  a 
plane  that  has  a  fixed  orientation;  (1v)  the  contact  normal  performs  In  the  pro¬ 
cess  of  meshing  a  parallel  motion;  (v)  the  motion  of  the  contact  ellipse  provides 
Improved  conditions  of  lubrication;  and  (vl)  the  gears  can  be  manufactured  by  use 
of  Gleason's  equipment. 
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